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Finite-size effects are investigated in tiie Gaussian model with isotropic and anisotropic short- 
range interactions in film geometry with nonperiodic boundary conditions (b.c.) above, at, and 
below the bulk critical temperature Tc. We have obtained exact results for the free energy and the 
Casimir force for antiperiodic, Neumann, Dirichlet, and Neumann-Dirichlet mixed b.c. in 1 < d < 4 
dimensions. For the Casimir force, finite-size scaling is found to be valid for all b.c. For the free 
energy, finite-size scaling is valid in 1 < d < 3 and 3 < d < 4 dimensions for antiperiodic, Neumann, 
and Dirichlet b.c, but logarithmic deviations from finite-size scaling exist in d = 3 dimensions for 
Neumann and Dirichlet b.c. This is explained in terms of the borderline dimension d* = 3, where the 
critical exponent l — a — f = (d— 3)/2 of the Gaussian surface energy density vanishes. For Neumann- 
Dirichlet b.c, finite-size scaling is strongly violated above Tc for 1 < d < 4 because of a cancelation 
of the leading scaling terms. For antiperiodic, Dirichlet, and Neumann-Dirichlet b.c, a finite film 
critical temperature rc,flim(i) < Tc exists at finite film thickness L. Our results include an exact 
description of the dimensional crossover between the d-dimensional finite-size critical behavior near 
bulk Tc and the (d— l)-dimensional critical behavior near TcfiimiL). This dimensional crossover is 
illustrated for the critical behavior of the specific heat. Particular attention is paid to an appropriate 
representation of the free energy in the region Tc,mni{L) < T < Tc- For 2 < d < 4, the Gaussian 
results are renormalized and reformulated as one-loop contributions of the </p* field theory at fixed 
dimension d and are then compared with the e = 4 — d expansion results at £ = 1 as well as with 
d = 3 Monte Carlo data. For d = 2, the Gaussian results for the Casimir force scaling function are 
compared with those for the Ising model with periodic, antiperiodic, and free b.c; unexpected exact 
relations are found between the Gaussian and Ising scaling functions. For both the d-dimensional 
Gaussian model and the two-dimensional Ising model it is shown that anisotropic couplings imply 
nonuniversal scaling functions of the Casimir force that depend explicitly on microscopic couplings. 
Our Gaussian results provide the basis for the investigation of finite-size effects of the mean spherical 
model in film geometry with nonperiodic b.c. above, at, and below the bulk critical temperature. 
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I. INTRODUCTION AND SUMMARY 

Critical phenomena in confined systems have remained 
an important topic of research over the past decades. 
Much interest has been devoted to systems confined to 
film geometry which are well accessible to accurate ex- 
periments, e.g., measurements of the critical specific heat 
and of the critical Casimir force in superfluid films near 
the A transition of ^He and ■^He-'^He mixtures [U |2] and 
in binary wetting films near the demixing critical point 
[3]. To some extent, these phenomena have been repro- 
duced by Monte Carlo (MC) simulations of lattice mod- 
els in finite-slab geometries [4H6]. While progress has 
been made in the theoretical understanding of these phe- 
nomena above and at the bulk critical temperature Tc of 
three-dimensional systems [7^^^ , there exists a substan- 
tial lack of knowledge in the analytic description of three- 
dimensional systems in film geometry below bulk T^, ex- 
cept for the case of periodic boundary conditions (b.c.) 



'Electronic address: bkastening@matgco.tu-dar mstadt.de| 
'^Electronic address: vdohm@physik.rwth-aachen.de 



|17j . except for the study of qualitative features of the 
critical Casimir force [IH] , and except for the study of dy- 
namic surface properties |19] . Also for two-dimensional 
systems in strip geometry, only a few analytical results 
have been known for the critical Casimir force [101 l20l - f23] 
in the past. Analytic expressions for the Casimir force 
scaling functions of the two-dimensional Ising model are 
known for free and fixed b.c. |23j and only since very 
recently for periodic and antiperiodic b.c. |24j . On the 
other hand, to the best of our knowledge, no complete 
analytic results for the free energy finite-size scaling func- 
tions are available for the elementary Gaussian model in 
strip and film geometries, respectively, in two and three 
dimensions for nonperiodic boundary conditions. 

There are several reasons for this lack of knowledge. 
One of the reasons is that realistic b.c, such as Dirich- 
let or Neumann b.c. for the order parameter, imply con- 
siderable technical difficulties in the analytic description 
of finite-size effects below bulk Tc even at the level of 
one-loop approximations. A second reason is the di- 
mensional crossover between finite-size effects near the 
i/iree-dimcnsional bulk transition at Tc and the two- 
dimensional film transition at the separate critical tem- 
perature TcfiimiL) < Tc of the film of finite thickness L. 
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An appropriate description of this dimensional crossover 
constitutes an as yet unsolved problem even for the sim- 
plest case of film systems in the Ising universality class 
with unrealistic periodic b.c. A third reason is the inap- 
plicability of ordinary renormalized perturbation theory 
to the (fi^ model in two dimensions (i.e., either at fixed 
dimension d = 2 or within an e expansion in 4 — e dimen- 
sions extrapolated to e = 2) because of the large value of 
the fixed point of the renormalized four-point coupling at 
d = 2. No special reason exists, on the other hand, as to 
why no attention has been paid in the literature to the 
Gaussian model in d = 3 film or d = 2 strip geometries 
with several different boundary conditions, although this 
model is exactly solvable and does provide valuable and 
interesting information on various aspects of the free en- 
ergy and the Casimir force, as we shall demonstrate in 
this paper. A short summary of our main results is given 
below. 

(i) Gaussian model as the basis for the mean spheri- 
cal model : The exactly solvable mean spherical model 
(MSM) [25] has played an important role in the analy- 
sis of finite-size effects near critical points where, how- 
ever, the free energy and the critical Casimir force have 
been studied, for a long time, only for periodic b.c. 
A calculation of the critical Casimir force in the MSM 
for nonperiodic b.c. was performed recently |27| . with a 
few results in film geometry in 2 < d < 3 dimensions. 
Clearly these results need to be extended to a more com- 
plete investigation. A serious shortcoming of the MSM 
is the pathological behavior of the surface and finite-size 
properties in d > 3 dimensions [25, 28J with logarith- 
mic deviations from scaling in d = 3 dimensions. Such 
logarithms were also found in the Casimir force and the 
free energy [571 US] . A profound understanding of these 
pathologies is important for the appropriate interpreta- 
tion of the deviations from finite-size scaling in the MSM. 
It was suggested earlier [30] that the pathologies in the 
MSM should be attributed to the effective long-range in- 
teraction induced by the constraint. The earlier analyses 
for nonperiodic b.c. (see |26j). however, were restricted 
to integer dimensions d = 3, 4, . . .. A more recent study 
[31j of the full continuous range of 2 < d < 4 dimensions 
revealed the absence of pathologies for d < 3 and identi- 
fied the origin of the nonscaling features for d > 3 as a 
consequence of the properties of the ordinary Gaussian 
model with short-range interactions. The crucial point 
is that the MSM can be considered as a Gaussian model 
with a constraint and that there exists a borderline di- 
mension d* = 3 in the Gaussian model above which the 
Gaussian surface energy density has a nonuniversal finite 
cusp at bulk Tc. This cusp causes all nonscaling effects 
for d > 3, while for d* = 3 the logarithmic divergence 
of the Gaussian surface energy density explains the loga- 
rithmic deviations from scaling in the three-dimensional 
MSM [M] . Both pathologies enter the MSM through the 
Gaussian surface terms of the constraint equation. The 
long-range interaction induced by the constraint does not 
yet introduce a nonuniversal parameter but it is rather 



the combination with the borderline dimension d* = 3 
of the Gaussian model with short-range interactions that 
is the origin of the nonuniversal nonscaling features for 
d > 3. The analysis of [31] was restricted to the regime 
t > with t = {T - Tc)/Tc for Dirichlet b.c, without 
considering the critical Casimir force. Our goal is to fully 
explore the finite-size critical behavior of the free energy 
and the critical Casimir force of the MSM both above and 
below Tc for five different b.c. and to properly explain 
the expected deviations from finite-size scaling in three 
dimensions as well as to study the scaling functions for all 
b.c. in 2 < d < 3 dimensions. It is our conviction that this 
goal must be based on a profound analysis of the Gaus- 
sian model as a first step, before turning to the MSM. The 
appropriateness of this strategy was demonstrated earlier 
in [31 . In the present paper we perform this first step. 
Our main results that will be relevant to our forthcom- 
ing analysis of the MSM are as follows, (a) Our results 
provide an exact description of the dimensional crossover 
from the d-dimensional finite-size critical behavior near 
bulk Tc to the (d— l)-dimensional critical behavior near 
7c, film: which is illustrated in Sec. |VIII for the critical 
behavior of the specific heat. This dimensional crossover 
will constitute the basis for describing the corresponding 
crossover from bulk Tc to T ^ for d < 3 in the MSM. 
(b) Our exact calculation includes nonnegligible logarith- 
mic non-scaling lattice effects in d = 3 dimensions for the 
case of Neumann b.c. and Dirichlet b.c. that have not 
been captured by the method of dimensional regulariza- 
tion used in Ref. [8j. Such effects will be important for 
the interpretation of the logarithmic nonscaling behav- 
ior in the d = 3 MSM model, (c) For the case of mixed 
Neumann-Dirichlet (ND) b.c, a strong power-law viola- 
tion of scaling is found in general dimensions 1 < d < 4 
that has an important impact on the scaling structure 
of the free energy density in a large part of the L~^^'^- 
t planes of both the Gaussian model and the MSM and 
that is expected to imply unusually large corrections to 
scaling in the (p"^ theory. 

(ii) Gaussian model scaling functions as one-loop 
renormalization- group (RG) scaling functions: There is 
another important reason for studying finite-size effects 
of the Gaussian model. After appropriate renormaliza- 
tion, the Gaussian results for the free energy, Casimir 
force, and specific heat can be reformulated as one-loop 
contributions of the (p* field theory. From previous work 
[52) it is known that, within the minimal subtraction 
scheme in d = 3 dimensions [33], the one- loop bulk am- 
plitude function of the specific heat provides a reasonable 
approximation above Tc and that the one-loop finite-size 
contributions for Dirichlet b.c. [TTl |T2] yield good agree- 
ment with specific- heat data p]l34| of confined ^He in film 
geometry above and at the superfluid transition. This 
suggests to determine the one-loop results for the free 
energy and the critical Casimir force within the minimal 
subtraction scheme at fixed dimension d and to compare 
these results with e = 4 — d expansion results at £ = 1 
[i H m [H], with recent MC data [H H and with 
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the recent result of an improved d = 3 perturbation the- 
ory [17] in an L| X L slab geometry with a finite aspect 
ratio p = LlL\\ — 1/4. As suggested by the earlier suc- 
cesses IT71 [55] , the minimally renormalized ip^ theory 
at fixed d is expected to constitute an important alter- 
native in the determination of the Casimir force scaling 
function in comparison to the earlier e expansion ap- 
proach [HI ini [H] ■ It is one of the central achievements 
of this paper that our d — i one-loop RG results shown 
in Fig. [5] below indeed support this expectation. 

(iii) Casimir force scaling functions in two dimensions: 
Most of our Gaussian results are valid in 1 < d < 4 
dimensions. This permits us to study the interesting 
case d — 2 and to compare it with the exact results 
of the two-dimensional Ising model [221 HSl IMJ- As a 
totally unexpected result we find (in Sec. VI) surpris- 



erties in d > dimensions, and give a short account of 



ing relations between the Casimir scaling functions of 
the Gaussian model with periodic (antiperiodic) b.c. and 
those of the Ising model with antiperiodic (periodic) b.c. 
Our comparison between these models also identifies the 
magnitude of non-Gaussian fluctuation effects in the two- 
dimensional Lp'^ model for several b.c. 

(iy) Nonuniversal anisotropy effects: It has often been 
stated in the earlier and recent literature [S] [TOl dH [HI 
12111261137] that the critical Casimir force scaling functions 
are universal, i.e., "independent of microscopic details". 
In view of these claims we briefly study the case of a sim- 
ple example of anisotropic couplings, i.e., two different 
nearest-neighbor couplings Jy and Jj_ in the horizontal 
and vertical directions, respectively. Our exact results 
for the Gaussian model show that these anisotropic cou- 
plings imply nonuniversal scaling functions of the Casimir 
force that depend explicitly on Jy and Jj_ for all b.c, as 
predicted by Chen and Dohm [351 ISSHiO] and recently 
confirmed by Dantchev and Griineberg f5P for the case 
of antiperiodic b.c. in the large-n limit for 2 < d < 4. 
In particular, we verify for all b.c. the exact relation 
[3S1I1I] Aaniso = (Jj./J||)^'^"^^/^Aiso between the Casimir 
amplitudes of the isotropic and anisotropic film system 
within the d-dimensional Gaussian model. We also ex- 
tend this kind of relation to the two-dimensional Ising 
model for periodic and antiperiodic b.c. in the form 
Aaniso = (?o,_l/Co,|| ) Aiso, where ^o,_l and ^o,|| are the 
correlation-length amplitudes perpendicular and parallel 
to the boundaries of the Ising strip. For the case of free 
b.c. at Tc, such a relation was found earlier by Indekeu et 
al. [20]. It would be interesting to test such nonuniver- 
sal anisotropy effects by MC simulations for the critical 
Casimir force, in addition to those for the critical Binder 
cumulant 42J. 

As a general remark we note that the Gaussian model 
does not have upper or lower critical dimensions; for this 
reason many of our results are valid for arbitrary d > 
except for certain integer d where logarithms appear (at 
even integer d for bulk properties and odd integer d for 
surface properties). 

The outline of our paper is as follows. In Sec. [llj we 
define our model, review the relevant bulk critical prop- 



what effects arise if the model is anisotropic. In Sec. Ill 
we consider the film critical behavior in 2 < d < A di- 
mensions. In Sec. |IV[ we derive and discuss the singular 
contributions to the free energy density in 1 < d < 4 di- 
mensions. In Sees. [VHVll the Casimir force is considered, 
in Sec. |VII| our results are reformulated as one-loop RG 
results of the (p"^ field theory and are compared to other 
RG and MC results, while in Sec. VIII we focus on the 
specific heat and its crossover from d to d — 1 dimensions. 
The Appendix is reserved for details of our calculations. 



II. GAUSSIAN MODEL IN FILM GEOMETRY 

A. Lattice Hamiltonian and basic definitions 

We start from the Gaussian lattice Hamiltonian (di- 
vided by k^T) 



n = d'^ 



ro 
2 



(2.1) 



with = Yla=ii^^'')'^ ^^"i '^i*'^ couplings Jx,x' 
between the continuous ri-component vector variables 

Sx — {Sx^ , • ■ • , Sx^^ ) on the lattice points a; of a d-dimen- 
sional simple-cubic lattice with lattice spacing d. The 
components vary in the range —00 < S'i"^ < +00. 
Unless stated otherwise, we shall assume an isotropic 
nearest-neighbor ferromagnetic coupling Jx.x' — J > 0, 
Jx,x' = for |£C — a;'| > a. In the discussion of our re- 
sults we shall also comment on the case of anisotropic 
short-range interactions Jx.x' with a positive definite 
anisotropy matrix A [36^ as defined in Eqs. (2.431 and 



(2.52) below. The only temperature dependence enters 
via To = agt = aoiT — Tc)/Tc, ao > 0, where Tc is the 
bulk critical temperature. We assume 
lattice points in a finite rectangular box of volume V — 
Ly"^ X L = Ma''', where Ly = TVyd and L = Nd are 
the lattice' extension in the d — 1 "horizontal" directions 
and in the one "vertical" direction, respectively. Thus 
we have N layers each of which has N^~^ fluctuating 

variables. The lattice points are labeled hy x = (y, z) 
with y = {yi, . . . ,yd-i)- We assume periodic b.c. in the 
horizontal (y) directions. As we shall take the film limit 
iVy — >■ cxD, the relevant b.c. are those in the vertical (z) 
direction. The top and bottom surfaces have the coordi- 
nates zi = d and z^ = L, respectively. It is convenient 
to formulate the vertical b.c. by adding two fictitious lay- 
ers with vertical coordinates z^ — Q and zjv+i = L -\- a 
below the bottom surface and above the top surface, re- 
spectively, for each value of the d — 1 horizontal coordi- 
nates. Then we may define periodic (p), antiperiodic (a), 
Neumann-Neumann (NN), Dirichlet-Dirichlet (DD), and 
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Neumann-Dirichlet (ND) b.c. by 



p 




^^N + l = 


Szl , 




(2.2a) 


a 










(2.2b) 


NN 


Szo 




•^Zjv + l 


— Szh J 


(2.2c) 


DD 


Szo 


= 0, 


•^Zjv + l 




(2. 2d) 


ND 


SzQ 






= 0, 


(2.2e) 



where we have omitted the y coordinates. We use the 
representation 



(z,(7)n4^(2/nft), (2.3) 



Ml ^(2, 9m) 



UL^(2,9m) 



(NN) , ^ 



(ND). ^ 



COS qrnZ = 1 

\pl COS (7„i z 

cos qmZ — COS ■ 
^/2 sm 9m -2' 

\/2 COS 9„j z 

cos 9m — COS ■ 

•\/2sin9„iZ 



1 J cos9,„(z - |) = 1 
\V2cos9m(z- 1) 



iV + 1 



sin9mZ 



— cos 9m (z - I) 



m = 

1 < TO< iV/2 

m = iV/2 

iV/2 < m < TV - 1 

< TO < (TV - l)/2 

m^[N - l)/2 

(iV- l)/2 < m < TV - 1 



9m = 



TO = 
TO = 1, 

TO = 0, . . 



.,iV-l 

,7V- 1 



TO = 0, 1 



27r(TO H 



■Km 

7r(TO + 1) 
L + a ' 

7r(TO + i) 



(2.4a) 



(2.4b) 

(2.4c) 

(2.4d) 
(2.4e) 



with the Fourier aniphtudes Sp^q and the complete set 

(t) 

Uj^ of real orthonormal functions, where, for the d — \ 

horizontal directions, the u^^\z, qm) are used with the 
replacements L — ^ Ly, z — > y^, and 9m — >■ Pi.m; — 
27rTOi/i||. The TO = N/2 mode for periodic b.c. (the 
TO = (A'^— 1)/2 mode for antiperiodic b.c.) is only present 
if N is even (if N is odd). The above mode functions 
are equivalent to those in [25l |3TJ |43], where complex 
mode functions for periodic and antiperiodic b.c. have 
been used instead of our real mode functions. 



nalized Hamiltonian 



The functions (2.4) satisfy the orthonormality condi- 
tions 



^ULiZj,qm)uLiZj,qm') = Sm,m', (2.5a) 

Zj 

Zj,qm)uL{zj',qm) - (2.5b) 



with Zj = ja, j = 1, . . . ,iV. For the case of isotropic 
nearest-neighbor couplings J > 0, this yields the diago- 



p,d-l 



p,q 



Jp.d-i = ^ (1 - cospiO.) 

4 = 1 

J„ = ^ (1 - cos 9a) . 



d-l 



(2.6) 



(2.7a) 



(2.7b) 



Equations (2.7 1 reflect the cubic anisotropy of the lat- 
tice. The lowest modes have p — and are homoge- 
neous (90 = 0) for periodic and NN b.c, whereas they 
are z-dependent with 90 = '^/{L + d) for DD b.c. and 
9o = 7r/(2L+a) for ND b.c. For antiperiodic b.c, there is 
a twofold degeneracy of the lowest modes with 90 = ir/L 



and 97V- 1 = — tt/L -I- 27r/a, since Jq^ = Jq 



This has 



important consequences for the behavior of the free en- 
ergy and specific heat near the film critical temperature, 
see Sees. |IIH |IV A[ and |VIII| below. A corresponding 
twofold degeneracy of the ground state is known for the 
mean spherical model with antiperiodic b.c. |41) . 
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We note that the boundary conditions assumed in ( 2.2 ) 
do not depend on any nonuniversal parameter. They are 
conceptually simple and represent only a small subset of 
a large class of more complicated boundary conditions. 
The latter may exist in the presence of an anisotropic 
lattice structure whose symmetry axes are not orthog- 
onal to the boundaries but have skew directions rela- 
tive to the boundaries. Such more complicated systems 
(which, however, belong to the same bulk universality 
class as standard spin models — such as Ising models with 
nearest-neighbor couplings on simple-cubic lattices) in- 
deed exist, e.g., among real magnetic materials with a 
non-orthorhombic lattice structure. Models of such sys- 
tems may also arise after a shear transformation has been 
performed to an isotropic system [36' if the original lat- 
tice model has non-cubic anisotropics. In this case the 
transformed boundary conditions depend on the original 
anisotropy parameters and therefore give rise to nonuni- 
versal finite-size effects. We shall come back to such skew 
nonuniversal boundary conditions in the context of the 



Sec. III). For antiperiodic, DD, and ND b.c, roc,fiim(i) 



discussion of two-scale factor universality in Sec. II C 
The dimensionless partition function is 



Z{t,^,L) = 


n 

.3/, 2 




exp( 


-H) 


Loo a(2-<^)"/2 




n 




exp( 








M 


27r 






52 (ro + JpM-i 







where we have used that, due to the orthonormality of the 
UL^ the linear transformation Sy^z ^ Sp^q has a Jacobian 

\{dSy,zldSp,q)[^ l. 

The film limit is defined for d > 1 by letting Ly — >■ 00 
while keeping L finite. In this limit the Gaussian free 
energy per component and per unit volume divided by 
k^T is given for t-q > roc,fiim(i) by 

/(t,L) = -i hm -^\nZ{t,^,L) 

1 1 ^ fid-i) 

" ^''^^^^ ^2L^J ^'^^ iro+Jp,d-i+Jq)] , 

(2.9) 

where jjf^^^ = HtiV-Jya ^^/(^Tr)- A /iZm critical 
point exists at ro = 7'0c,fiim(i), where the argument of 
the logarithm on the right hand side of (2.9) vanishes for 
p = and q — qq. 

As a shortcoming of the Gaussian model, the bulk crit- 
ical value roc = and the film critical value roc,fiim(-^) are 
independent of d and n, and no low-temperature phase 
exists. Furthermore, the Gaussian roc is not affected by 
lattice anisotropics, in contrast to rQcfiimiL) , which de- 
pends explicitly on the anisotropic couplings J^.x' (see 



is negative, thus the free energy (2.9) exists for negative 
values of vq in these cases. The region roc,fiim(i) < ^0 < 
will be of particular interest for the study of the mean 
spherical model below the bulk transition temperature 
|44j . The film critical behavior of the Gaussian model 
will be discussed in more detail in Sec. IIIII 

The bulk limit is obtained by letting L — > 00, 
L-'Eq ^ = I-^i~^dq/{2n). The bulk free energy 
density per component divided by k^T is, for t > 0, 

/b(i) = /(<,c^) 

1 1 f^'^) 

= H2TT) + In [a^ (ro + Jk,d)] ■ (2.10) 



In the long- wavelength limit, the cubic anisotropy does 
not matter and J^.d = 2Jk^ + 0{k'^) becomes isotropic 
which justifies to define a single second-moment bulk cor- 
relation length ^ above T^, 



lim 



1 E^.^'i^-^ns^s^') 



L-X^ 2d 
The latter is given by 
^=(2J/ro)l/2=^o^-^ 



(2.11) 



^o = (2J/ao)^ iy^l/2. 

(2.12) 



In the presence of NN or DD b.c, there are surface free 
energy densities per component 2f^^\t) and 2f^f\t) for 
i > as defined by 



/sf(t) 



1 

2l 



lim {L[f{t,L)-Mt)]}. 



(2.13) 



In the presence of ND b.c, the total surface free energy 
density per component is 



p(D) 



(2.14) 



For periodic and antiperiodic b.c. there exist no surface 
contributions. 

For small t > 0, the bulk and surface free energy den- 
sities will be decomposed into singular and nonsingular 
parts as 

/b(i) = /b,s(i) + /b,ns(t), (2.15a) 

/sf(i) =/sf,s(i)+/rf,ns(t), (2.15b) 

where /b,ns(i) and /sf,ns(0 have an expansion in positive 
integer powers oit. For small t and large L, it is expected 
EH HH ESj that, for the Gaussian model ([2T]) in film 
geometry, the free energy density can be decomposed as 



fit,L)^Ut,L) + U,{t,L), 



(2.16) 



Usit, L) = A,ns(i) + L-'[fl°Ut) + fZM. (2-17) 
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where "top" and "hot" refer to the top and bottom sur- 
faces of the film. In the absence of logarithmic bulk sin- 
gularities [21], i.e., for d ^ 2 and d ^ 4, and in the 
absence of logarithmic surface singularities |31) . i.e., for 
d ^ 3 (or periodic or antiperiodic b.c), the singular part 
is expected to have the finite-size scaling form [46] 



(2.18) 



with a nonuniversal parameter Ci. For given b.c, the 
scaling function J-{x) is expected to be universal only 
within the subclass of isotropic systems but nonuniver- 
sal for the subclass of anisotropic systems of noncubic 
symmetry within the same universality class |36l I38| , see 
( 2.44 )-( 2.58) below. A convenient choice of the scaling 



varia 



3le X is 



X = t{L/^^f/\ 



(2.19) 



i.e., Ci = ^^'^ . The bulk singular part 

hAi)^YdC\ d>0, 2,4,6,..., (2.20) 



see Sec. |IIB| belo w, wi th a universal bulk amplitude Yd is 
included in Eq. (2.181 through = limj_j.oo a^"''" J^(i) 



for 1 < d < 4, d 7^ 2. For the surface free energy density, 
(|2.18|) implies 



l-d 



(2.21) 



a universal surface 



amplitude 



with 
limj_ 

For fixed i > and large L it is expected [H] [21 [Ml iZ] 

that the free energy density can be represented as 



f{t,L)^Mt)+L-'[fJ^{t) + f^r{t)] 



L-^Gix) 



0(e 



(2.22) 



In (2.22), is the exponential hulk correlation length in 



the direction of one of the cubic axes [SHI HZ] 



e..(^arsinhA 
\~a 2^ 



(2.23) 



Its deviation from ^ for finite d causes scaling to be vio- 
lated [3S1[17] for fixed t > and large L > 24^^/5."^, i.e., 
X > 576(^/d)'*. We recall that this scaling violation is a 
general consequence of the exponential structure of the 
excess free energy for large L at fixed ^ and is a lattice 
(or cutoff) effect that is predicted to occur not only in 
the Gaussian model but quite generally in the (p^ lattice 
(or field) theory for systems with short-range interactions 
[36[ |47] . This effect is different in structure from addi- 
tional nonscaling effects that occur in the presence of 
subleading long-range (van der Waals type) interactions 
[361148]. 

In the absence of long-range interactions, no contribu- 
tions ~ L~™ with m > 1, m d should exist in (2.22) 



for film geometry. The representation (2.22) separates 



the finite-size part ~ L~'^ from the surface parts ~ L~^. 
The latter do not contribute to the Casimir force scaling 
funct ion X {x) to be d is cussed in Sec. [V] 

If (|2.18|) and (|2.20|)-(|2.22|) are vahd, the connection 



between J^, Agt, and Q is, for i > 0, 

J^ix) = Ydi'''' + (A*°P + A^[')i^''-^> + g{i). (2.24) 

In Sec. jlVj w e shal l examine the range of validity of the 
structure of ( |2.18[ ), ( |2.22p , and ( [2.24^ for the Gaussian 
model for various b.c. and calculate the scaling functions. 



In Sec. VIII we shall also discuss the specific heat (heat 



capacity per unit volume) divided by /cb 

C{t,L) =dU{t,L)/dT, 



(2.25) 



where U{t,L) — ~T^df(t,L)/dT is the energy density 
(internal energy per unit volume) divided by k^, with 
the singular bulk part 



f/b,s(i) - -T.Co'^'diyYdC^'-'^^^''- 



(2.26) 



The surface part of the energy density is Usi{t) = 
—T^dfsf{t) /dT, with the singular part 

UsfAt) = -Tc^o'^'id - l)i.A,fr (2.27) 



In ( 2.26[ ) and (2.271 we have used the hyperscaling rela- 
tion dif — 2 — a, with the Gaussian exponent 



a 



(4-d)/2, 



d < 4, 



In the presence of NN, ND, and DD b.c. 
deviations from the scaling structure of (2.18|, (|2.21[), 



(2.28) 
logarithmic 



(2.24), and (2.27) are expected for the Gaussian model 
in the borderline dimension d* = 3 [31] because of the 
vanishing of the critical exponent 



1-a-iy^ (d-3)/2 



(2.29) 



of the singular part of the surface energy density ( |2.27[ ). 
(This is similar to the logarithmic deviations for systems 
with periodic b.c. 45^ at d = 4, where the specific-heat 
exponent a vanishes.) In this case, J^{x) and ^sf do not 
exist, but G{x) and X{ x) rem ain well defined. The posi- 
tivity of the exponent ( [2.29 1 for d > 3 implies a nonuni- 
versal cusp that is responsible for the nonscaling features 
in the MSM for d > 3 [3I]- 

Moreover, logarithmic deviations from the structure of 
(2.20) and (2.26) are expected for the Gaussian model in 



the borderline dimension d = 2 because of the vanishing 
of the critical exponent 



l-a = (d-2)/2 



(2.30) 



of the singular part of the bulk energy density (2.26) 



B. Bulk critical properties 

In contrast to real systems with short-range interac- 
tions, the Gaussian model has a bulk phase transition at 
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ro = for any dimension d > including d = 1. In the 
following we present both the singular and nonsingular 
parts of the bulk critical behavior of the free energy since 
they will be needed in the context of the mean spherical 
model in a subsequent part of the present work 44 . The 
exact result for the bulk free energy density for rp > in 
d > dimensions is 



with the singular part 



/b,s(i) 



ln(C/a) 



d = 2. 



(2.40) 



1 

25? 



where fo = r^d? /{2J) and 



In 



J 



dy 

y 



(2.31) 



(2.32) 



with 



B{y) EE e-yioiy), 



(2.33) 



and where Iq is a Bessel function of order zero, Io(z 



dimensions is derived as 

n-d/2) 



TT ^ dip exp(z cos (fi). From the large-y behavior ( B12 ) 
of B{y), the universal amplitude of (2.201 in d > 0, d ^ 
2,4,6 



The logarithmic structure is related to the vanishing of 
1 — a for d = 2, see (2.30). In Sec. VI we shall compare 
the Casimir force scaling function of the Gaussian model 
with that of the two-dimensional Ising model. This com- 
parison will be restricted to the regime T > Tc- Cor- 
respondingly, we comment here on the Ising bulk free 
energy only for this case. For T > Tc the bulk correla- 
tion length of the d = 2 Ising model is, asymptotically, 
^ = ^0+^ with V = 1. In terms of this length, the 
singular part of the bulk free energy density of the d=2 
Ising model (on a square lattice with lattice spacing a) 
has the same form as given by (2.401 but with a nega- 
tive amplitude — l/(47r) instead of l/(47r) for the d = 2 
Gaussian model. 



In contrast to the universal power-law structure (2.20) 
for d > 0, c? ^ 2,4,6, . . ., the logarithmic structure (2.40) 



2(47r)'i/2 ' 



(2.34) 



with Ya = — (127r) ^. The nonsingular bulk part /b,ns(i) 
has an expansion in integer powers of fo, 

1 



contains the nonuniversal microscopic reference length 
a. Other reference lengths are expected for other lattice 
structures, whereas the amplitude l/(47r) is expected to 
be universal. The choice of the amplitude of such ref- 
erence lengths is not unique but in our case the lattice 
spacing d appears to be most natural for the cubic lat- 
tice structure. (Due to the artifact of the Gaussian model 



/b,„s(t) = /b.„s(0) + ^ [/i7~o + 0{fl)\ , (2.35) and the d = 2 Ising model that ^"^ ~ t and ^"^ ~ t 
where 

/b,n.(0) = /b(0) = \n-+WM , (2.36) 



respectively, are analytic functions of a different choice 
cd with c 7^ 1 as a reference length would yield a different 
decomposition into singular and nonsingular parts.) 



/i = <! 27o 



dy \B{yf ~ (2^y)-'^/2 



, < d < 2, 
d>2. 



(2.37) 



C. Isotropic and anisotropic continuum 
Hamiltonian 

For the purpose of a comparison with the results of tp'^ 
field theory we shall also consider the continuum version 



with the generalized Watson function 



1 



Wa{z)^W'^{z)^- dye-^y/^B{yY. (2.38) 



of the Gaussian lattice model (2.1) for an n-component 
vector field <p{x). For the choice 2 J = 1 the isotropic yi^ 
Hamiltonian reads 



In order to appropriately interpret the critical be- 
havior of the three-dimensional system in film geome- 
try in Sees. |III||VIII below it is important to first con- 
sider the hulk critical behavior in two dimensions. While 
1^2(0) = 4G/7r with Catalan's constant G w 0.915966 
is finite, both Yd = -l/[47r(d - 2)] -t- 0((d - 2)°) and 
Wdi'S) = l/[27r(d - 2)] + 0{{d - 2)°) diverge as d -> 2+. 
However, the sum of the respective contributions to the 
singular and the nonsingular part of the free energy re- 
mains finite and we obtain the bulk free energy per unit 
area 



d^x 



fQ 2 
2 ^ 



i V — 

2^ \dx^ 



2\2 



(2.41) 



with some cutoff A in k space. The field ip{x) = (p{y, z) 
satisfies the various b.c. that are the continuum analogues 
8J of Eqs. (2.2). In Sec. VII our Gaussian results based 
on (2.1), in the limit a — > will be renormalized and 



reformulated as one-loop contributions of the minimally 
renormalized ip* field theory at fixed dimension 2 < d < 4 
El [33] based on Hfidd, (|2.41|), in the hmit A 00. The 



/b(i) = /b(0) 



ln(^/d) ln2-l 



471^2 



167rJ 



(2.39) 



role played by the d = 3 RG approach will be to change 
the Gaussian critical exponent = 1/2 to the exact crit- 
ical exponent v at d = 3 entering the correlation length 
^ which appears in the scaling argument of the scaling 
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functions of the renormalized ip^ field theory. This will 
then justify to compare the resulting one-loop finite-size 
scaling functions of the Casimir force in d = 3 dimensions 
with MC data for the three-dimensional Ising model [B], 
with higher-loop e = 4 — d expansion results at e = 1 
[H [15] , and with the recent result of an improved c? = 3 
perturbation theory |17) in an L| x L slab geometry with 
a finite aspect ratio p = L/L\\ = 1/4. 

We shall also consider the anisotropic extension of 



(2.41) 



^ field 
aniso 



= / A 



Agfi dip dip 
2 dxa dxR 



+uo{ipy 



(2.42) 



The expression for the symmetric anisotropy matrix A = 
{Aap) in terms of the microscopic couplings Jx,x' of the 
lattice Hamiltonian (2.1), is given by the second mo- 
ments I36ll39l 



'^/3 ) Jx,x' 



(2.43) 



In the case of isotropic nearest-neighbor couplings J on 
a simple-cubic lattice we have simply Aq,^ — 2J5ap- In 
general, Aap is non-diagonal and contains d{d+ l)/2 in- 
dependent nonuniversal matrix elements. 

The relation between the finite-size critical behavior of 
isotropic and anisotropic systems was recently discussed 
in detail for the case of a finite rectangular geometry with 
periodic b.c. [MlEHlIll]- It was shown that the relation 
between the anisotropic and isotropic critical behavior is 
brought about by a shear transformation. In real space, 
this transformation is described by the matrix product 
^-i/2[/^ with an orthogonal matrix U that diagonalizes 
A according to A = UAU^^, where A is a diagonal ma- 
trix whose diagonal elements are the eigenvalues of A. 
This transformation causes a nonuniversal distortion of 
the rectangular shape to a parallelepipedal shape, of the 
simple-cubic lattice structure to a triclinic lattice struc- 
ture, and of the periodic b.c. along the rectangular sym- 
metry axes to periodic b.c. along the corresponding skew 
lattice axes of the triclinic lattice. The general structure 
of the scaling form of the free energy density is expressed 
in terms of the characteristic length L' = V'^l'^ where V 
is the finite volume of the parallelepiped (see Eqs. (1.3) 
and (4.1) of |36]). This is, however, not directly appli- 
cable to our present model with film geometry with an 
infinite volume and with various b.c. Furthermore, a 
significant difference occurs in film geometry due to the 
existence of a film transition temperature that is affected 
by anisotropy for the cases of antiperiodic, DD, and ND 
boundary conditions. Thus anisotropy effects in film ge- 
ometry deserve a separate discussion. In particular, we 
shall compare our results with those of Indekeu et al. 
|20j . who studied an anisotropic Ising model on a two- 
dimensional infinite strip. 



In [5S] it was found that, for 2 < d < 4 in the large-n 
limit of the theory above bulk in film geometry with 
periodic b.c, the universal structure ( 2.18| ) is replaced by 



,{t,L;A) = L-''[{A-'),,]-''/^:FUiL/e)'^n, 

(2.44) 



where J-iso is the scaling function of a film system de- 
scribed by the isotropic ip"^ theory with ordinary peri- 
odic b.c, but where the scaling argument contains the 
transformed length 



L=[{A-'),d]'/'L, 



(2.45) 



and where i^' is the b ulk correlation length of the isotropic 



system. In (|2.44|), {A ^)dd denotes the dth diago- 



nal element of the inverse of the reduced matrix A — 
A/{detAy/'^. In [3F the simplicity of the structure of 
(2.44) was attributed to the large-n limit. In general one 
expects that /s,aniso(ij ^) is expressed in terms of the 
scaling function of an isotropic system that has trans- 
formed boundary conditions which are not identical with 
those of the original anisotropic system. For a brief dis- 
cussion of such boundary conditions see the paragraph 
before Eq. ([ZS]) in Sec. [lIA) 

A simplifying feature of film geometry is that the shear 
transformation preserves the film geometry except that 
the original thickness L is transformed to a different 
thickness L. In general, the length L appearing in the 
scaling argument of J^so in (2.44) is not the transformed 



thickness L but rather the distance between those points 
on the opposite surfaces in the transformed film system 
that are connected via the periodicity requirement |38| : 
this distance is measured along the corresponding skew 
lattice axis. The correctness of this geometric interpreta- 
tion can be seen as follows. Let Xd = zhe the unit vector 
in the ^-direction, i.e., orthogonal to the film boundaries. 
Then L is the length of the vector L obtained by trans- 
forming the vector Lx^i, i.e., 

L ^ X-^^^ULxd, (2.46) 

and therefore 

L = \\-'^^ULxd\ = \xlU-^\-^'^\-^'^Uxd\^'''L 
= \xlU-^X-^Ux,\"^L = \xlA'^xa\''^L 
= [{A-^)mY'^L, (2.47) 



in agreement with (2.45). A corresponding statement 



holds for antiperiodic b.c. 

As we show in Appendix |a] the thickness L of the 
transformed isotropic film is given by 



L = (detA-Vdct[[A-i]])i/2 



(2.48) 



where the (d— 1) x (d— 1) matrix [[A^-'^]] is obtained by 
removing the dth row and column from A^^. 

It is possible to express (2.44) in terms of the single 
length L by rewriting 

/,,a„iso(i, L- A) = (det A)-i/Vs,iso(i, L), (2.49) 
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(2.50) 



Thus, apart from the geometric factor (det A)^^/^ that 
describes the change of the volume of the primitive cell 
under the shear transformation, /s^aniso is given, in the 
large-71 limit, by the free energy fs,iso{t, L) of an isotropic 
film with an effective thickness L ^ L with ordinary pe- 
riodic b.c. We conjecture that the structure of (2.49) 



with (2.50) is exactly valid also for the Gaussian model 



with periodic b.c. A similar structure is expected to be 
valid for the Gaussian model with antiperiodic b.c. ex- 
cept that the scaling argument should be expressed in 
terms of t rather than ^ in order to capture the regime 
Tcfiim < T < Tc^buik- Furthermore, the effect of the 
anisotropy on Tc^fiim needs to be taken into account (see 
Sec, [m] below). 

A nontrivial situation exists in the case of NN and ND 
b.c. because Neumann b.c. involve a restriction on the 
spatial derivative perpendicular to the boundary which, 
after the transformation, turns into a derivative in a 
skew direction not necessarily perpendicular to the trans- 
formed boundary. Thus the isotropic film system still 
carries the nonuniversal anisotropy information of the 
original system both in its changed thickness and in 
the nonuniversal orientation of its transformed boundary 
conditions. Thus both nonuniversality and anisotropy 
are still present at the boundaries of the transformed sys- 
tem. The same assertion applies to periodic and antiperi- 
odic b.c. Clearly, since boundary conditions dominate 
the finite-size critical behavior at Tc where the correlation 
lengths extend over the entire thickness of the film sys- 
tem, the above reasoning implies that universality is not 
restored by the shear transformation in spite of internal 
isotropy (in the long-wavelength limit) of the transformed 
system away from the boundaries. In other words, even 
this internal isotropy of a confined system does not en- 
sure the universality of its critical finite-size properties 
because of the nonuniversality contained in the bound- 
ary conditions. In the light of these facts we consider as 
incorrect the recent assertion by Diehl and Chamati [50] 
that "the critical properties of an anisotropic system can 
be expressed in terms of the universal properties of the 
conventional (i.e., isotropic) theory". 

More specifically, even after the shear transformation, 
the finite-size effects of the transformed isotropic system 
still depend, in general, on ci(d-|-l)/2-|-l nonuniversal pa- 
rameters (see Eqs. (1.3)-(1.5) of [3S]), contrary to the 
hypothesis of two-scale factor universality [22l |4^ . This 
multiparameter universality is fully compatible with the 
general framework of the RG theory [36]. Technically, 
these parameters enter through the transformed wave 
vectors k' of the isotropic system, thus the dependence 
of finite-size properties on A^j^ cannot be eliminated by 
the shear transformation as demonstrated explicitly for 
the example of periodic b.c. in Eq. (2.22) of [36]. We con- 
clude that there is no basis for complying with the tra- 
ditional picture of two-scale factor universality according 
to the suggestion "to define two-scale factor universality 



only after the transformation to the primed variables (of 
the isotropic system) has been made" [50]. This sugges- 
tion would be applicable only to hulk properties of the 
transformed system. 

A special case is the case of DD b.c. (vanishing order- 
parameter field if at the boundaries) or free b.c. (no 
condition on the fluctuating variables at the boundaries) 
since these b.c. are invariant under the shear transforma- 
tion and therefore do not violate isotropy. In particular, 
these b.c. do not introduce any nonuniversal parameter. 
Nevertheless, even in this case there is a nontrivial shift 
of Tcfiim of the film critical point of systems in the or- 
dinary {d,n) universality classes (for n = 1, d > 2, for 
n ^ 2, d > and for n > 2, d > 3 ) due to anisotropy. 
For the special case d = 2, n = 1, however, i.e., for a 
system of the Ising universality class on an infinite strip 
of finite width, there is no separate "film" transition and 
thus no analog of a finite Tc^m-m > exists. This con- 
ceptually simplest case was studied by Indekeu et al. [50] 
as will be further discussed below. One may conjecture 
that for DD b.c. the structure of (2.44) is valid also for 
the d-dimensional Gaussian model where, however, the 



length L in (2.44) is to be replaced by L, (2.48) 



An open question remains as to what extent the struc- 
ture of ( 12.44[ ) with \2Ab\ (and correspondingly of ( 12.58[ ) 
below) is valid even for the full y'* model with finite n in 
d dimensions and even for real film systems. It would be 
interesting to explore this problem theoretically as well as 
by means of MC simulations for a variety of anisotropic 
spin models in film geometry with various b.c. and vari- 
ous anisotropics. 

The situation becomes particularly simple if the ma- 
trix A is diagonal in which case the original simple-cubic 
lattice of the anisotropic system is distorted only to an 
orthorhombic lattice of the isotropic system that still has 

-1/2 r 



dd 



a rectangular structure. Then we have L = L = 
see App. |XJ In the following we confine ourselves to this 
simple case. 

We consider only two different nearest-neighbor inter- 
actions J|| and Jl in the "horizontal" and "vertical" di- 
rections. This corresponds to replacing Eqs. (2.7) by 



4J| 



Jp.d-1 = V! (1 ^ cospid) , (2.51a) 

i=l 

4:J± 

J„ = -3— (1 — cosga) , (2.51b) 

in which case A is given in three dimensions by 

J|| \ 

A = 2 I J|| . (2.52) 
J_L / 

In this case we must distinguish two different correlation 
lengths and ^j^. For the Gaussian model they are given 
by 

Cll =eo,||i"", ^0,11 = (2J||/ao)", ^^ = 1/2, (2.53a) 

a = ^o.±r^ co,± = (2J±/ao)^ ^. = 1/2. (2.53b) 
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The existence of two different correlation lengths implies 
the absence of two-scale factor universality [36l [38l |39] . 
As a consequence, all bulk relations involving correlation 
lengths have to be modified and all finite-size scal- 

ing functions are predicted (38i to become nonuniversal 
as they depend exp licitly on the ratio Jl/J|[. 
For the example (2.52), we obtain 



A 



(J±/J||)i/3 
{JJJ\\V^' 

{J±/Jii)-'/' 



(2.54) 

for d = 3 and iA-^)dd = {J±/ J\\Y^~''^''' and {A-^)dd = 
(2Jl)^^ for general d. For the isotropic Gaussian model, 

— 1/2 

we have simply ^' — rQ_ (compare Eq. (B16) of [55]). 



Then the scaling form (2.44) becomes 



ys,aniso i.^: i J\\ t J-L ) 

= L-''iJjJlO^'-'^^^T,.o{tiL/^o,±)'h- (2.55) 

Since 



(2.56) 



according to (2.531, this relation can be written as 



(2.57) 

where the finite-size scaling function J-aniso of the 
anisotropic system, considered as a function of the single 
scaling variable t{L/S^Q^±_)^^^ , is nonuniversal, 

^a„iso(i(i/eO,±)'/"; J||,./±) 

= (eo,±/eo,||)'"'-^iso(t(Veoa)'/^), (2.58) 

as it depends on the nonuniversal ratio Ja_/J\\ through 
the factor (^o,±/Co.|| )'^^^- (For the d = 2 Ising model (see 
(2.60) and Sec. VI B), this factor depends on J± and 
J|| separately.) As a consequence, also other thermody- 
namic quantities have a corresponding finite-size scaling 
structure. This was recently confirmed for the case of 
antiperiodic b.c. in the large-n limit in 2 < c? < 4 di- 
mensions [H]. So far no explicit verification of (2.55) 



has been given for systems with surface contributions. 
In App. [b] we shall verify that (2.551 holds within the 
Gaussian model for all b.c. in 2 < d < 4 dimensions, in- 
cluding those involving surface terms, in the temperature 
range where finite-size scaling holds. The consequences 
for the Casimir force scaling functions will be discussed 
in Sec.|v] Eq. (2.55) is not directly valid for the free en- 
ergy density in d = 2 dimension since the bulk part has 
a logarithmic structure, see (2.40), but we have verified 



that it is valid for the excess free energy density and for 
the Casimir force scaling form of the d = 2 anisotropic 
Gaussian model for all b.c. (see Sees, [v] and VI) 



The issue of nonuniversality of finite-size amp 



itudes of 



studied earlier in the work by Indekeu et al. PD]. In this 
paper an anisotropic Ising model on an infinitely long 
two-dimensional strip with free b.c. in the vertical direc- 
tion was considered. This corresponds to our geometry 
for the special case d ~ 2 with DD b.c. As noted above, 
this is a particularly simple case as no distortions of the 
b.c. arise even if the anisotropic couplings correspond to 
a nondiagonal anisotropy matrix. Furthermore, there ex- 
ists no analog to a "film" transition at finite width of the 
infinite strip below the two-dimensional "bulk" critical 
temperature T^. since there exists no singularity in an 
effectively one-dimensional system with short-range in- 
teractions. For the present case of interest, i.e., for the 
case of two different nearest-neighbor couplings in the 
horizontal and vertical directions, the Ising Hamiltonian 
(divided by ksT) of Indekeu et al. [20J contains ferro- 
magnetic nearest-neighbor couplings denoted by Ki and 
K2 which in our notation correspond to 2/3 J|| and 2/3 Jj_, 
respectively, with a lattice spacing d — 1. 

The authors derived an exact relation between the am- 
plitudes Aaniso and Aiso of the free energies at criticality 
of the anisotropic and isotropic Ising strips of the form 



sinh(4/3,J^) 
sinh(4/3,J,i) 



1/2 



Since 



[sinh(4/3eJ±)/sinh(4/3,J||)]i/2 ^ 



(2.59) 



(2.60) 



is the ratio of the amplitudes of the correlation lengths 
perpendicular and parallel to the Ising strip [20) . 



Eq. (2.591 can be written as 



(2.61) 



This is the same structure as given in (2.58) for d ~ 2. 

It was also shown that the ratio Co,_l/^o,|| can be inter- 
preted as a geometrical factor that arises in a transforma- 
tion of lengths such that isotropy is restored ^\ . This is 
in complete agreement with the analysis presented here 
and in Refs. 1361 and p39 . Nevertheless, in spite of the ex- 



act relation (2.591, it is clear that restoring isotropy does 
not imply "restoring universality" [20] since the finite-size 
amplitude Aaniso of the original anisotropic lattice model 
depends explicitly on the microscopic couplings J|| and 
J± ■ 

We note that the dependence of Co,±/Co,|| on the 
nearest-neighbor couplings Jy and Jj_ is a nonuniversal 
property that has a different form in d = 2 dimensions 



the free energy with respect to coupling anisotropy was 



for the Gaussian model on the one hand (see (2.56)) and 
for the Ising model on the other hand (see (2.60)) ■ The 
latter is not captured by heuristic arguments based on a 
mapping of a d = 2 lattice spin model on a continuum 
model as seen from Eq. (6.5) of Ref. [50] . 



III. FILM CRITICAL BEHAVIOR 

In the following we briefiy discuss the film critical be- 
havior of the d-dimensional Gaussian model which we 



11 



need to refer to in Sec. IIVI Here we confine ourselves to 
2 < d < 4. 

First we consider the isotropic case. For finite i, the 
film critical point is determined by tq = 7'oc,fiim(^) with 
?'Oc,fiim(-Zj) — for periodic and NN b.c, whereas 



roc,fiim(i) = -(4J/52)[1 - cos(goa)] < 0, 



(3.1) 



with qo = TT / L for antipcriodic b.c, qq — 7r/(i + a) for 
DD b.c, and go = 7r/(2L + a) for ND b.c, respectively. 
For large L/d, roc,fiim(-^) = — 2j7r^/i^ for antiperiodic 
and DD b.c. and roc,fiim(^) = -2j7r'^/{AL'^) for ND b.c. 
Correspondingly, the film critical lines are described, for 
large L, by 

t,,mmiL) = [Te,mm(i) - Tel/Tc = -TT^ (^0/^) l/"- (3.2) 

for antiperiodic and DD b.c. and by 

tc,mm{L) - -(^/2)2(eo/L)i/'^ (3.3) 

for ND b.c, in agreement with finite-size scaling. For the 
shape of the film critical lines see Figs. [2] and [3] below. 

Near Tc.mmiL) there exist long-range correlations par- 
allel to the boundaries. A corresponding second-moment 
correlation length ^fiim(?'Oj L) may be defined by 



y,z,y',z'iy y')^{Sy,zSy',z' 



2(rf-l) Ey,z,y',z'{Sy.zSy> 



(3.4) 



(The summation over all z, z' corresponds to a kind of 
averaging over all horizontal layers.) 
Define a length 



2J 



1/2 



(3.5) 



For periodic and NN b.c, where tqc, film (i) — 0, we ob- 



tain just as in the bulk case (2.12) the exact relationship 

1/2 

Cfilm = i = £, ^ (2J/ro) , which is independent of L. 
For antiperiodic, DD, and ND b.c. and arbitrary L/d, we 
obtain ^gim — I only in the limit where L. 

At finite L, the free energy per unit area divided by 
k-aT is defined as 



fmn.{ro,L)^Lf{t,L). 



(3.6) 



We expect that for £ ^ L (which is equivalent to the 
condition ^sim S> L mentioned in [8]), the film critical 
behavior corresponds to that of a bulk system in d — 
1 dimensions. Taking into account (2.201, this would 



imply that the singular part fmm.s has the temperature 
dependence for 2 < d < 4, d ^ 3, 



/film,s(i) — Yd-lC 



d-l^film 



d^3. 



(3.7) 



where the dimensionless universal amplitude Y^-i is de- 
fined by (2.20). We indeed confirm this expectation for 



all b.c. except for antiperiodic b.c. whose lowest mode 



has a two-fold degeneracy as noted already in Sec. II A 
above. This causes a factor of 2 in the corresponding 
relation 



d ^ 3, antiperiodic b.c. 

(3.8) 



For d = 3, the expected structure of /fiim.s is less obvi- 
ous because of the logarithmic dependence of the corre- 
sponding bulk quantity (2.39) in d = 2 dimensions. For 
£ ^ L we obtain 

/film(?'0, L) = /film(?'Oc,film, L) + /film.s + 0(^*0 " ?'Oc,film), 

(3.9) 

where it is necessary to specify the singular part /fiim.s 
separately for the various b.c. 



Ap) _ 

•J film,s 


^efiiLln(&hn/i), 


(3.10a) 


/•(a) _ 
J film,s 


^^filmln(efilm/i), 


(3.10b) 


.(NN) _ 
•f film,s 


^^filmln(&lm/VZa), 


(3.10c) 


„(DD) _ 
film,s 


^efiiLin(&.mavvi^/^), 


(3.10d) 


„(ND) _ 
J filin,s 


^^filmln(&lm/i), 


(S.lOe) 



as wifi be derived in Secs. |IVAl |IVC| and |IVD| below 



Both microscopic and macroscopic reference lengths 
may appear in the logarithmic arguments depending on 
the b.c. (Our decomposition is such that no logarithmic 
dependencies on a or L appear in the nonsingular part of 
/film proportional to vq — roc,fiim-) By contrast, the am- 
plitude l/(47r) appears to have a universal character, in 
agreement with ( |2.40 ), except for the factor of 2 for an- 
tiperiodic b.c. The expressions for /fiim(''Oc,fiim, L) at the 
critical line T = rc_fiim(i) are nonuniversal and depend 
on the b.c. 

For the anisotropic case, we consider only two dif- 
ferent nearest-neighbor interactions Jn and J± as de- 



scribed by (2.511 and corresponding different bulk cor- 
relation lengths 1^11 and as described by (2.53). In the 



paragraph containing Eqs. (3.1)-(3.3) the only necessary 



changes are a replacement of J by J± and of by Co._l- 

Define the film correlation length, now called ^fii„ 
by (3.4) and lengths £± and £ii by 



ilm,amso ; 



2J, 



1/2 



(3.11) 



For periodic and NN b.c, we obtain, in close analogy 
to the isotropic case, the exact relationship ^fiim.aniso — 

1 /2 

^11 ~ ^11 ~ (2*^11 /''o) 1 which is again independent of L. 
For antiperiodic, DD, and ND b.c. and arbitrary L/d, we 
obtain ^fiim.aniso = ^|| Only in the limit where i± 3> L. 

The considerations of the paragraphs containing 
Eqs. (3.6)-(3.10) translate one to one to the anisotropic 
case if ^gim is replaced by ^fiim,aniso and the condition 
^ ^ L is replaced hy £± ^ L. 
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IV. FREE ENERGY IN 1 < d < 4 DIMENSIONS 



A. Periodic and antiperiodic b.c. 



In the following we present exact results for the asymp- 
totic structure of the finite-size critical behavior of the 
Gaussian free energy density near the bulk transition 
temperature for large L/a 1 in the isotropic case. 
These results include both the bulk critical behavior 



( |2.20[ )-( |2.40[ ) for L oo at fixed t > and the film criti- 
cal behavior (p7t-( |3.10| ) for T Tcfiin,{L) at fixed finite 
L. Thus our results provide an exact description of the 
dimensional crossover from the d-dimensional finite-size 
critical behavior near bulk Tc to the (d— l)-dimensional 
critical behavior near Tcfiim of (the isotropic subclass of) 
the Gaussian universality class. Our scaling functions 
T are analytic at bulk Tc for antiperiodic, DD, and ND 
b.c, in agreement with the general discussion given in 
Sec. VII of Ref. 8^ . Our Gaussian results go beyond the 
corresponding one-loop results of Ref. [5j in the following 
respects: (i) Our exact calculation includes nonnegligi- 
ble logarithmic non-scaling lattice effects in d = 3 di- 
mensions for the case of NN and DD b.c, whereas these 
effects are not captured by the method of dimensional 
regularization used in Ref. [5]. (ii) For the case of ND 
b.c, a strong power-law violation of scaling is found in 
general dimensions 1 < d < 4 that has an important 
impact on the scaling structure of the free energy den- 
sity in a large part of the L~^/'^-t plane of the Gaussian 
model and that is expected to imply unusually large cor- 
rections to scaling in the theory, (iii) Our represen- 
tation of the scaling functions is directly applicable to 
the region Tc,fara{L) < T < Tc for antiperiodic, DD and 
ND b.c, whereas the representation of Ref. [8] is applica- 
ble only to r > Tc, apart from a few results in Sec. VII 
of Ref. JS]. (iv) We study the approach to the critical 
behavior near Tc^fiim and compare it with the critical be- 
havior of a (d— l)-dimensional bulk system; this compar- 
ison confirms the unexpected factor of two of the leading 
universal amplitude for the case of antiperiodic b.c. that 



was presented in our Sec. Ill as a consequence of the 
twofold degeneracy of the lowest mode, (v) Our analy- 
sis includes, for all b.c, the exponential nonscaling part 
of the excess free energy due to the lattice-dependent 
nonuniversal exponential bulk correlation length (2.23) 



that was not taken into account in [8]. (vi) Our analysis 
also includes the d = 2 scaling functions of the finite- 
size part of the free energy that provide the basis for the 
Casimir force scaling functions in d = 2 dimensions to be 
discussed in Sec. [v] and VI (it is only the d = 2 bulk part 
of the free energy that exhibits a logarithmic deviation 
~ ln(^/d) from scaling, see Sec. II B I; the case d = 2 was 
not discussed in [8]. 

Because of the special role played by the borderline di- 
mension d* = 3 for the surface properties of the Gaussian 
model it is necessary to distinguish the cases without sur- 
face contributions (periodic and antiperiodic b.c.) from 
those with surface contributions (NN, DD, and ND b.c). 



For periodic and antiperiodic b.c, the finite-size scal- 
ing structure of ( |2.18[ ) and ( |2.22 ) is confirmed. For d 7^ 2 
we find (see Appendix [b|) the finite-size scaling functions 

J-(P)(i) = 4p^(5) + Ydi'^'^, i > 0, (4.1a) 



with the universal bulk amplitude from (2.34), where, for 
2/>0, 



1 /"^ 



(4.2a) 



27r 



'^^[K{z/4)-K{z)]-^n/z 



(4.2b) 



Eq. (4.2b) is valid for 2 < d < 4, while for 1 < d < 2 



the subtraction of the term ^ y/irz inside the curly brack- 
ets has to be omitted. The function K{z) is defined by 
K{z) = X^n^oo exp(— n^z), which converges rapidly for 
large z. It may be expressed in terms of the third ellip- 
tic theta function i?3(u, e~^) [51] via K{z) = ?93(0, e~^). 
It satisfies the relation K{z) = yj -k j zK{t{^ j z^ with the 
expansion K{z) = \/t^/ z'^^^_^ exp(— ri^7r^/z), which 
converges rapidly for small z. 

The function J"(^)(i) is regular at i = in agreement 
with general analyticity requirements, whereas J^'^^\x) is 
nonanalytic at a; = due to the film critical point. 



Eqs. (4.1) include the singular parts of both the 
bulk critical behavior (x — >■ 00) and the film criti- 
cal behavior {x — -> for periodic b.c. and 
X + Tr'^ = i^/Cfiim ~^ for antiperiodic b.c). The latter 
is obtained from the singular parts of the small-y expan- 
sions for y > 

l(P)(y) =l(P)(0) + r,_iy('^-i)/2 + O(y,//2), (4.3a) 
I^;\y)^I^;\0) + 2r,_iy('^-i)/2 + 0{y,z//'h (4.3b) 

for d 7^ 3, while for d = 3 

(4.4b) 

Contrary to the naive expectation based on universality, 
the amplitudes of the leading singular and y\ny 

terms of (4.3) and (4.4), respectively, differ by a factor 
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of two for periodic and antiperiodic b.c. as already men- 



tioned in Sec. Ill These terms yield the right hand sides 



of ^Jj, pl^, ( |3.10a| and ( pUOb ). 

Comparison of (2.22) and ( |4.ip ieads to the finite-size 
parts for i > 



g(^)(i)-2i-'^g(p)(45)~g(p)(i), 



(4.5a) 
(4.5b) 



where (4.5b) follows from (B6a). Eqs. (4.5) remain valid 
for d = 2. For d = 3, Eg. ( |4.5b[ ) agrees with Eqs. (9.3) 
for iV = 1 of Ref. [8 . 

The representation of our results differs from that of 
[H], where Eqs. (6.8) provide an integral representation 
of Q^^\x) and 0^^'>{x). Both representations have the 
same expansions in terms of modified Bessel functions, 
see Appendix [C| which suggests that, for i > 0, in- 



and g(-^{i) = eWp„,(y+) 



deed g^P^x) = e^^l,iy, 
with the identification y^^- = \fi. Our representation of 

T^'^\x), e^^P^i), and Q^^'^iyx) in terms of J^p^ has the 
advantage that it is directly applicable to the bulk crit- 
ical point at i = 0, whereas the integral representation 



(1) 

+por 



(?/+) and e^^^p„,(i/+) given in Eqs. (6.8) of [S] 



of e 

require an extra small-2/+ treatment of the divergent in- 
tegrals so that after multiplication with the prefactor yi^ 
finite results are obtained. More importantly, the repre- 
sentation of T^^^ [x) in terms of X^f' has the advantage 
that it is valid also for i < including the film critical 
point at X — — tt^, whereas the integral representation 

of 6^aper(2/-i-) in Eqs. (6.8) of [5^ is not suitable for an 
analytic continuation to the region i < 0. 

A representation of J-^^^ {x) valid for all i, > — tt^ may 
also be extracted from the result (3.26) with (3.27) in 
Ref. fU] for the singular part of the excess free energy 
of the mean spherical model with antiperiodic b.c. in 
film geometry. After omitting the term oc Xt, restor- 
ing the bulk contribution by removing the term cx yj^^^, 
and replacing y ^ x + tt'^, the last two terms within the 
curly brackets of Eq. (3.26) in Ref. l4] may be sho wn to 
be eq uivalent to the integral representation (4.1b) with 
( |4.2b| of J'('')(i). 

The universal finite-size amplitudes at Tc are 

J-(P)(0) = g<^P\0) = -Tr-'^/^T{d/2)C{d), (4.6a) 
^(a)(o) = g(a)(o) = (1 - 2i-'^)^-'^/2r(rf/2)C(d), (4.6b) 
which agree with the corresponding iV = 1 amplitudes 



ApCT and Aaper, respectively, in Eq. (5.7) of 8_ (up to a 
sign misprint there for periodic b.c). 

At fixed i > the results for ^(p) and g^""^ yield the 
large-i approach to the bulk critical behavior 



.(1) 



1 



X > 1, (4.7) 



where the upper (lower) sign refers to periodic (antiperi- 
odic) b.c, see the paragraph around (B13). For suffi- 
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FIG. 1: Scaling function J"'''' (x), ( 4.8b I, of the free energy of 



the Gaussian model in three dimensions with antiperiodic b.c. 
for X > —tt'^ (solid line). For i > are also shown the bulk 
part Ysx^^'^ (dotted) and finite-size part Q^^\x) (dashed). 



scaling form (4.7) must be replaced by an exponential 



nonscaling form |36j which is obtained from (4.7) by re- 
placing the exponential argument — -s/i by — i/^o, where 

is the exponential correlation length (2.23). 

In d = 3 dimensions the scaling functions (4.1 ) can be 
expressed as 



^(P)(i)=g(P)(J)--Lj3/2^ 
iZTT 

J-(-)(i) =^(-)(£)- ^^3/2 



127r 



i > 0, 



X > — TT 



with the finite-size parts 
1 

'2n 
1 

2^ 



g'-pHs:)=- 



Li3(e"^") + ViLi2(e-^= 
Li3(-e-^) + \/£Li2(-e-^") 



(4.8a) 
(4.8b) 

(4.9a) 
(4.9b) 



where Li„(z) are polylogarithmic functions (see Ap- 
pendix Id]). With the identification y-^- ^ x wc find 
that g'^ix) and g^'^\x) agree with Q+pcr{y+) and 
Q+apei:{y+) in Eqs. (9.3) of |H], respectively, but that 
the representation of 0+apci (y+) is more elaborate than 
that of C/ ('')(£). It is understood that in (4.8b) for 
— TT^ < i < 0, th e fun ction g'-'^^x) means the analytic 
continuation of (4.9b) to i < which is complex; to- 



i^/2/(127r), however, the 



gether with the complex term 
right-hand side of ( 4.8b ) becomes real and analytic for all 
X > —TT^ with a finite real value 



J-(^)(-7r2) = - 



03) 
2tt 



-0.191313, 



(4.10) 



ciently large L at fixed t > 0, however, the exponential 



see Appendix IeI For x > 0, the d — 3 scaling functions 
g'--^\x) and g'^{x) will be shown in SecMtogether with 
the corresponding scaling functions X^^^'Jx) and A"'^^)(i) 
of the Casimir force. 

In Fig. [l] we show the scaling function J"(^)(x), (4.8b), 
of the Gaussian free energy density in three dimensions 
for antiperiodic b.c. including the range for negative x 
down to the film transition at i = — tt^. It would be in- 
teresting to compare this result with the corresponding e 
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expansion result at e = 1 which, however, is not available 
in the literature so far. 



B. NN and DD b.c. in d 7^: 3 dimensions 

For NN and DD b.c. there exist well-defined surface 
free energy densities for i > in d > 1 dimensions. They 
are given by (see Appendix [b|) 



dy 



Biy)"-' [e-'y - 1] e 



-Vfo/2 



1 



1 



y 

VKd-i(ro)-W^rf-i(4 + fo) 



r^B{yr''[e-'y + l-2Biy) 

Jo y 



(4.11a) 



-yro/2 



8a' 



d-l 



-Wd_i(fo) - Wd-i{4. + fo) + 2Wd{fo)\ , 

(4.11b) 



with fo defined after Eq. (2.31). The result for f^f'^ 
agrees with /surface in Eq. (67) of [3T]. For d 7^ 3 and 
small t > 0, the singular parts are 



A. 



(D) 



r(-i) & 



(4.12a) 

+ 0(r''^+'^), (4.12b) 



4(47^)^/2 

with the universal surface amplitudes 

^If' = = J^'d-i. l<rf<5, d^3, (4.13) 

in agreement with Eq. (6.3) in [B] and the remark about 
the surface contribution in the last paragraph on page 
1910 of [H] as well as with Eqs. (76) and (88) in [3T]. The 
nonsingular parts are for r = N, D 



/i2s(^) = ftp (0) - + 0(r~g), (4.14) 



a 

with the nonuniversal constants 

/•OO 



1 

16 7o 



dy 



B{y) 



-l) + (2^2/)(i-'^)/2 



>0, 
(4.15a) 



1 1'°° 

bf^ ^-j^ dy[B[yY-'[e-^y + 1 - 2B[y)] 

- (27ry)(i-'')/2 + 2(27ry)-'^/2} > 0. (4.15b) 



for b^^"* are connected by analytical continuation in d. 



Eq. ( |4.15b[ ) holds for 1 < d < 2, while for 2 < d < 3 (where 
&^^^<0) the addition of 2(2TTy)-'^/'^ inside the curly 
brackets has to be omitted. For 3 < d < 5 (where Vf'^ > 
again), additionally the subtraction of (27ry)(^~'^)/^ has 

to be omitted. These integral expressions for V^J^'' are 
connected by analytical continuation in d, as already 
noted in [3T], where closely related integral representa- 
tions of 6^°^ were given in Eqs. (77) and (89), valid for 
2 < d < 3 and 3 < d < 5, respectively. For d -^• 3, both the 
amplitudes and the coefficients Up diverge, while 
the nonuniversal co nstant s f^'' (0) < and f'^'' (0) > re- 
main finite, see Sec. 



IV C 



below. For d ^ 2, Uf'^ diverges 



but the corresponding term in (4.14) combines with the 



subleading term in (4.12b) to give a finite contribution 
to fT' 



The terms in (4.12) agree with the corresponding 
contributions in Eq. (6.3) and Appendix C o f [8). The 
sum of (|4.12b| and (|4.14|) for r = D b.c. with (14. 131), and 



V'J'^ as given in Ref. 31 agrees with Eqs. (75)-(77) and 



(87)-(91) of [3T]- In ( |4.12b[ ) we have included a s ingula r 
term of order Such a term does not exist in ( |4.12a |. 
Although this term is subleading compared to the leading 
singular ^^^'^ term, it becomes a leading singular term 



for ND b.c. (to be discus sed in Sec. IV D below), where 
the terms of (4.12a| and (4.12b) cancel because of 
(pSl). 



For NN and DD b.c. the finite-size scaling structure 
of (|2.18[) and \2.22\ is confirmed for d 7^ 3. We find the 



finite-size scaling functions (see Appendix [BJ) 

-f 2A([)(i-f c7r2)('^-i)/2, i>-cn\ (4.16) 

with c = for r = N and c = 1 for r = D, with the 



universal bulk amplitude Yd from (2.34), and where 

(4.17a) 



Tr\y) = 2-^I^p{Ayl 



Q (y) = 

r(p) 



dz (^/z)(^+l)/2x 



— \/TTZ + z| , 

(4.17b) 



withl^^^ from (4.2a). Eq. (4.17b) is valid for 3<d<4, 



cal behavior (3.7) 



r2 = iVel 



Eq. (4.15a) holds for 1 < d < 3, while for 3 < d < 5 (where 
&|(^^<0) the addition of (27ry)(i-'')/2 inside the curly 
brackets has to be omitted. These integral expressions tained from the surface terms of Eqs. (4.16) and from 



while for 2 < d < 3 (1 < d < 2), the addition of z (the ad- 
dition of z and the subtraction of \/t^) inside the curly 
brackets has to be omitted. The function F^^^\x) is 
regular at i = in agreement with general analyticity 
requirements, whereas !F^^^^ (x) is nonanalytic at i; = 
due to the film critical point. 

Eqs. (4.16) include the singular parts of both the bulk 
critical behav ior (2.20) (i -> 00) and the film criti- 

-> for NN b.c. and 
The latter is ob- 



filn 



■> for DD b.c). 
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singular parts of the small-y expansions for y > 0, 

d^3, (4.18a) 
d^2,3. (4.18b) 



At fixed t > the results for g'^^^^ and yield 
the same large-L approach to the bulk critical behavior 



fit,L)-Mt) = 



L'^ 2(47r)('i-i)/2 

i > 1, (4.22) 



We note that, according to (4.13), the surface amplitudes 
Asf of the d-dimensional film system have the same d- 
dependence as the bulk amplitude Yd-i of the (d— 1)- 
dimensional bulk system, apart from a constant factor of 
±1/4. This imphes 



see the paragraph around (B13). Eq. (4.22) is in agree- 



ment with the result Eq. (72]Tn [31] for free (DD) b.c. 
For sufficiently large L at fixed t > 0, the exponential 
part of the scaling form (4.22) must be replaced by an 
expon ential nonscaling form j36] which is obtained from 
(4.22) by replacing the exponential argument —2\fi by 



2A 



(N) 
sf 



(D) 
sf 



(4.19) 



—2L/^c, where is the exponential correlation length 
( 2.23 ). The same remark applies to the exponential parts 



which explains how the y'^'^ ^^/^ terms on the right hand 
sides of (4.18) and the terms in (4.16) involving the sur- 



contained in the scaling functions that are presented in 
Eqs. (|425l), (li^il), and (l438| below. 



face amplitudes (4.13) lead to identical amplitudes 1^- 



for the film free energy in (3.7) for both NN and DD b.c. 



in agreement with the expectation based o n univ ersality. 

For the finite-size contribution ^ L^'^ in (2.22) we find 
the scaling functions for i > 



NN and DD b.c. in d = 3 dimensions 



For NN and DD b.c. in d = 3 dimensions, the vanish- 



= ir'(-)=2-^4P)(45), 



ing of the critical exponent (2.29) of the surface energy 



(4.20) 



density [31 causes logarithmic deviations ~ ln(^/a) from 
the sc aling structure of ( |2.18[ ). From ( |4.12[ ), ( |4.14[ ), and 
(4.15), we obtain for d — > 3 the singular and nonsingular 



where g^^^\i) agrees with Eq. (71) in |31j with 
X = \fi > 0. 

The representation of our results differs from that of 
[5], where Eqs. (6.8) and (6.6) provide an integral rep- 
resentation of ^('^^^ and Q''^^\ respectively. Both rep- 
resentations have the same expansions in terms of mod- 
ified Bessel functions, see Appendix [C] which suggests 
that, for i > 0, indeed g'^^^'>{i) = e+sji^sjiiy+) and 
g(DD)^~-j _ 8^Q o(2/+)i with the identification ?/+ — \fi. 
Our representation of J"(nn)(j)^ g'^^^')(x), and g(°°)(i) 
in terms of 2^^^^ has the advantage that it is directly ap- 



parts of the surface free energy density for small t > as 



/if^(t) = /if^(0) + 



ln(C/5) _ ^ 
167r^2 " 



+ o(f2,r'ine), 

+ O(f2,r'ln0, 
with the nonuniversal constants 



1 a 

24^1^^ 



5(D) 



(4.23a) 
(4.23b) 



plicable to the bulk critical point at i = 0, whereas the = lim {h^^^ 

integral representation of 0+sBSB(y+) ^^"^ ®+o 0(2^-1-) '^'^^ 



A 



(N) 



given in Eqs. (6.8) and (6.6), respectively, of [S] require 
an extra small-?/+ treatment of the divergent integrals 
so that after multiplication with the prefactor yi^ finite 
results are obtained. More importantly, the representa- 
tion of T'^^^ (i) in terms of I^^^ has the advantage 
that it is valid also for a; < including the film critical 
point at 2; — — TT^, whereas the integral representation 
of 6^Qo(y+) in Eq. (6.6) of [5] is not suitable for an 
analytic continuation to the region a; < 0. 
The universal finite-size amplitudes at Tc are 

J-(NN)(0) ^ g(NN)(o) ^ ^(DD)(0) = g(DD)(Q) 

= -(4^)-'^/2r(d/2)C(d), (4.21) 



16 
1 



dy B{yf [e-^y - l] 



-y/2 



27ry 



W^2(4) 



5 In 2 + 1 



327r 

&(°)^lim iUF^-A'^ 



-0.027653, 



1 

' 32^ 
(4.24a) 



^^j^ dy\B{yf [e-^y + l-2B{y)\ 



1 - e-y/^ 

2Try 



+ 



32n 



^^[W,{4)-2WM] + ^-^^^ 



-0.00199279. 

(4.24b) 



which agree with the corresponding iV = 1 amplitudes 
^SB,SB for NN b.c. and a'^^^^ for DD b.c. in Eqs. (5.7) 
and (5.6) of [5], respectively. 



The limit d — > 3 in (4.24) is independent of whether 



it is taken as d — >• 3_ or d — 3+, in agreement with 
Eqs. (79) and (92) of [3T] for the case of Dirichlet b.c. 
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The structure of the leading singular terms of (4.23a) and 



(4.23b) agrees with that of the two-dimensional result 
(2.39) but the amplitudes are different. For the same 



reason as in ( 4.12b ), we h ave inc luded t he sub leading 
term in (4.23b 1. Eq. (|4.23b|) with (|4.24b|) agrees 



with Eqs. (80)-(82) of [31 but here we give a simplified 
expression of b'-^'' as compared to Eqs. (81) and (82) in 



The singular surface contributions ( 4.23a[ ) and (4.23b) 
appear also in the resulting singular parts of the free 
energy densities for < > (see Appendix [b]) , 



L3 ■' 



(4.25) 



with "+" for T = NN and " for r = DD, and where 

g(NN)(-)^^(DD)(~)^2-3g(P)(4i) 



1 

16^ 



Li3(e-'^) + 2ViLi2(e-^^^) 



(4.26) 



Eq. ( |4.25| ) for DD b.c. agrees with Eq. (86) in [31] (there 
is a sign misprint in Eq. (85) of '3T]). With j/+ = 
Eq. ( |4.26| agrees with Eqs. (9.3) for TV = 1 in [8j. 



Because of the dependence of ln(^/a) on the nonuni- 
versal lattice spacing a, no finite-size scaling functions of 
fs^^\t,L) and fs^^\t,L) can be defined. The nonuni- 
versal surface terms ^ L^^ constitute the leading de- 
viations from the bulk critical behavior for large L at 
fixed t > 0. One may define "nonscaling regions" in the 
{S^q/ LY/^-t planes (see Figs. [2|a) and (b)) by requiring 
that these logarithmic terms are comparable to or larger 
than the scaling terms L^^Q{x). These nonscaling re- 
gions depend on a/^o ^^nd are shown for the example 
^/■Co = 1 as the shaded regions in Figs.[2];a) and (b). 

The logarithmic deviations from scaling are not present 
right at bulk Tc, where /s(0,L) = ^"^^(O) with the 
universal critical amphtudes ^(^^^(0) = ^(^^'(O) = 
— C(3)/(167r), in agreement with Eq. (9.2) for = 1 in 

For the remaining part of the discussion we need to 
distinguish the cases of NN and DD b.c. For NN b.c, the 
function fs^^\t,L) is not regular at ^ — a s it includes 
the film critical behavior (3.9) with (3.10c) for t — >■ 
at fixed L. To derive this behavior we use the small- a; 
expansion for i > 



T(NN)^~; 



(i) 



C(3) i(lni+21n2-l) i^/^ 
167r 167r ^ 12^ 



which implies 



(NN), 



L3 



X ln{xa/L) 
8^ 



(4.27) 



0{x) . 
(4.28) 



The second term yields (3.10c) because of i = i^/Cfii„ 



for NN b.c. 
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FIG. 2: Asymptotic part of the {^o/L)'-^''-t plane of the 
Gaussian model in three dimensions with isotropic short- 
range interaction in film geometry with NN (a) and DD (b) 
b.c. The solid lines indicate the film critical temperatures 
Tc.min{L) at finite L: (a) vertical line at t = for NN b.c, (b) 
Eq. ( |3.2[ ) for DD b.c. No low-temperature phases exist for 
T < Tcfiiin{L). Finite-size scaling is valid between the film 
critical lines and the shaded areas. The shaded areas (as de- 
fined in the text) are nonscaling regions that depend on fi/^o. 
Their shapes are shown here for the example fi/^o ~ 1. These 
shapes start at the origin with infinite slope. The film criti- 
cal lines in (a) and (b) have the same form as for the cases of 
periodic and antiperiodic b.c, respectively. In these cases sur- 
face terms are absent and only very small nonscaling regions 
exist due to the nonscaling exponential parts ~ exp(— L/^o) 
mentioned after Eq. (4.7 1. 



By contrast, the film critical point for DD b.c. is lo- 
cated at tc,fiim < 0, see (3.2), thus no singularity exists 
at f = for finite L for DD b.c, which implies that 

/i^^'' [t, L) should be regular at i = 0. This is indeed the 
case as shown in the f ollowing. The first three terms of 
fs^^\t,L) from (4.25) can be rewritten as 



1 

Z3 



t(dd).~ 



(i)- 



12^ 



:lna; x\n{L/a) 



167r 



87r 



(4.29) 



where now the logarithmic deviation from scaling appears 
in the form of ln(L/a) but the temperature dependence 
through i ~ i is regular at t = since 



(DD)/~ 



[x) 



x?!"^ ilni 



127r 



167r 



(4.30) 



is regula r at x = (see Appendix [e| . It is understood 



that in (4.30) for -tt^ < i < 0, the function 5(°°'(i) 
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means the analytic continuation to a; < as given by 
(4.26), which is complex; together with the complex 
terms —x 3/7(127r) + iln5/(167r), however, ( |4.30[ ) be- 
comes real and analytic for x > 



read for t > 



value at i = 



■ with a finite real 
see Appendix |E| The representation 
(4.29 ) has the advantage that it is valid down to i = — vr^ 
corresponding to the film critical point. For x — >■ — vr^ it 



includes the film critical behavior (3.91 with (S.lOd). To 
derive this behavior we use an expansion around x — — tt^ 
for X > - 



t(dd).~. 



(5) 



' (see Appendix |E| 

x^/^ xlnx 



127r Wtt 

C(3) + 27r2ln7r (5 + 7r2)[21n(i + vr^) - 41n7r - 3] 

167r 16n 

+ 0((i + 7r2)2), (4.31) 

which implies 



/i°°)(t,L) = /i°°)(ic,fihn,i) 

^ f {i + TT^)\n[{i + Tr^)L/a)] 



+ 0ix + 7T^) 



(4.32) 



The second term yields (3.10d) because of 



i + n^ ^ L^/^li^^ for DD b.c. 

For i > 0, the d — 3 scaling tunct 
g(°°)(£) will be shown in Sec.jvltogc 
responding scaling functions X*^^^ ^ 
the Casimir force. 



functions ^^'^^'^^(a:) and 
;ether with the cor- 
(i) and X(°°)(i) of 



D. ND b.c. in 1 < d < 4 dimensions 

For ND b.c, the leading terms of the singular parts 
of the surface free energies, i.e., the 0{£,^~'^) terms in 
( 4.12a[) and (4.12b) and the logarithmic terms in (4.23a) 
and (4.23b), cancel. Then the leading term of the singu- 



lar part of the total surface free energy density for t > 0, 



r(-i) a 



0(C 



-id+l)) 



(4.33) 



does not have the universal scaling form (2.21), but de- 



pends explicitly on a. The cancelation of the leading sur- 
face terms ^ ^^-d fQj. ]\fj) ^.^ .^^j^g already noted in Ap- 
pendix C of [H, w here , however , the next-to- leading sur- 
face term ~ (4.33), was not taken into account. In 



contrast to the weak logarithmic deviations from scaling 
in d = 3 dimensions according to (4.23), Eq. (4.331 con- 



stitutes a strong power-law violation of scaling (within 
the Gaussian model) that has an important impact on 
the scaling structure of the free energy density in a large 
part of the L~^^^-t plane. The resulting singular and 
nonsingular parts of the free energy density for ND b.c. 



/(fl^)(i,i) = /b,„.(t) 



L 



n-i) a 

4(47r)<i/2 



g(ND)(3i) 

(4.34) 



1 
L 



where + Uf'^ < and, for x > 0, 

e(ND)(^)^2-'^g(-)(45), 



(4.35) 



(4.36) 



with Q^'^^ from (4.5b) (see Appendix |B[). This result re- 

^(D)^ _ 



mains valid for d 



3, where lim^^jjl^^-' 



&(N) _^ j(D) ^i^j^ ^(N) and b^P) given by (4.24). For d ^ 2, 



the divergent terms of ( 4.34 ) and ( 4.35 1 combine to give a 
finite result, while ^ {x) remains finite and continues 
to provide the scaling function of the finite-size contribu- 
tion to the free energy. 

The representation of our results differs from that of 
Ref. [5], where Eqs. (6.8) provide an integral representa- 
tion of Cj(^P). Both representations have the same ex- 
pansions in terms of modified Bessel functions, see Ap- 
pendix [cj which suggests that, for i > 0, Q^^^\x) — 
®+o SB(y+)i "with the identification y+ ^ Our rep- 
resentation of Q^^^\x) in terms of Q^^\x) and thus x'f'^ 
has the advantage that it is directly applicable to the bulk 
critical point at i = 0, whereas the integral representa- 
tion of 9*^Q sB^y+) given in Eqs. (6.8) of fS' requires an 
extra small-y+ treatment of the divergent integral so that 
after multiplication with the prefactor y'\_ a finite result 
is obtained. More importantly, the related representa- 



tion of J^(^P) (i) in terms of provided in Eq. ( 4.40 ) 



below has the advantage that it is valid also for x < 
including the film critical point at a; = — (7r/2)^, whereas 

the integral representation of 0+osB(y+) ^<\s. (6.8) 
of [8j is not suitable for an analytic continuation to the 
regio n x < 0. 

In (4.34) the nonscaling structure of the surface term 
^ destroys the finite-size scaling form of fs^^\t, L) 
above in the regime where the surface term is compa- 
rable to or larger than the finite-size term L~'^Q'-^^^ {x), 
i.e., in the regime 



T \ d-i ~ 
L \ ^ > 



4(4^)^/5 



n-d/2) 



(4.37) 



with X = (L/^)'^/'' for i > (this is valid only for d ^ 2; 
for d = 2 the logarithm in the singular part of the bulk 
free energy causes deviations from scaling in any case). 
For d = 3 this regime is indicated by the shaded area 
in Fig. [3j This violation of finite-size scaling is signifi- 
cantly more important than that due to the exponential 
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FIG. 3: Asymptotic part of the {^o/Ly^'^-t plane of the 
Gaussian model in three dimensions with isotropic short- 
range interaction in film geometry with ND boundary con- 
ditions. The solid line indicates the film critical tempera- 
ture Tc,flim(i) at finite L according to Eq. (3.3 1. No low- 
temperature phase exists for T < Tcfiim- Finite-size scaling is 
valid between the film critical line and the shaded area. The 
shaded area is a nonscaling region that depends on fi/^o- Its 
shape (shown here for the example a/^o = 1) is determined 
by Eq. (4.371. The shape of the shaded area starts at the 



origin with zero slope. 



correlation length ^c, ( 2.23 1, which happens only for con- 
siderably larger L > 24^'^/a^. At fixed i > 0, the result 



(4.34) yields the large-L approach to the bulk critical 



behavior 

/f°)(i,L)-/b,3W = L-^ 

j(rf-l)/4 



r(-i) ~a 



2(47r)('^-i)/- 



-e-2^/5, i>l, (4.38) 



see the paragraph around (B13). Eq. (4.38) implies that 



for ^/a ^ 1 the estimate (4.37) for the nonscaling region 
in Fig. [3] can be replaced by L > ^^ln(^/a). 

The cancelation of the leading surface scaling terms in 
the Gaussian model does not persist in the tp"^ theory at 
0{u*) (two-loop order) as can be seen from Eq. (Dll) of 
[H]. Two-loop terms, however, are typically smaller than 
one-loop terms, therefore it is expected that the two-loop 
contributions to the scaling part are less important than 
ordinary 1-loop scaling contributions. This means that 
now corrections to scaling are expected to become con- 
siderably more important compared to the scaling part. 
This would imply a shrinking of the asymptotic region for 
the case of ND b.c. Thus we expect that a more careful 
analysis of future experiments or of MC simulations is 
required for systems with ND b.c. because of unusually 
large corrections to scaling. 

Above the shaded area in Fig. [3] the nonscaling sur- 
face term is negligible and the leading L dependence of 
fs^^\t, L) is, for d ^ 2, described by the scaling form 

/fD)(t,L) = L-'^^(ND)(j)^ i>-(7r/2)2, (4.39) 



_^(ND)(~) ^ j(ND)(^ ^ ^^/2^2) ^ Y,[S: + (^/2)2]'^/2 

^ 2-'ij-(a)(-4~) (4 40) 



and 



ir^(,)^2-'^4^)(4y), 



(4.41) 



with T^^'^ and X^''-' from (4.1b) and (4.2b), respectively. 
This r esult includes the film critical behavior (3.7) and 
(3.10e) for x — >■ — (7r/2)^ at fixed finite L. To derive this 



behavior we use an expansion around x — —{n/2) , 

d^3, (4.42) 

while for d = 3 



m^vMMll^^Oiff^). (4.43) 



167r 



The second terms on the right hand sides of (4.42) and 



(4.43), respectively, yield (3.7) and (3.10c) because of 



The universal finite-size amplitude at Tc is 

J-(ND)(o) = g(ND)(o) = (1 - 2i-'*)(47r)-'*/2p(^/2)C(d), 

(4.44) 

which agrees with the corresponding = 1 amplitude 
^O.SB in Eq. (5.7) of 



For d = 3 we combine (|4.8b|) and (|4.40|) to obtain 

(4.45) 



1 

12^" 



7(ND)/~N _ 



1 

16n 



Li3(-e-2^) + 2\/iLi2(- 



(4.46) 



With the identification = y/i we find that 

g(ND)j-~-j agrees with the more elaborate representation 
of 0+o,SB(2/-f^) provid ed by Eq. ( 9.3) in Ref. Because 
of the relations (4.36) andJ 4.40 ), t he si tuation is similar 
to that explained after ji^^md (jij]), thus J'(^°)(i) 
is real for x > — (7r/2)^ and an analytic function for 
X > — (7r/2)^, even though the analytic continuation of 
g(ND)j^~-j .j-p negative x becomes complex. For i > 0, 
the d = 3 scaling bmction fj('^^)(a;) will be shown in 
Sec. |V] together with the corresponding scaling function 
X^-^^{x) of the Casimir force. 

In the region where finite-size scaling is valid (see 
Fig. [3]) there exists a scaling function J^^^^\x) of the 
free energy density for x > —{n/2)^. Due to Eqs. (4.36) 



and ( 4.40 ) , a plot of this function in three dimensions can 
be obtained from the solid curve in Fig. [l] with appropri- 
ately rescaled axes (the same holds for Q^^^{x) and the 
bulk part Yai^/^). 
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V. CASIMIR FORCE 

The excess free energy density per component divided 
by fceT is defined by 



/ex(t,i)=/(t,i)-/b(i), 



(5.1) 



where fhit), ( 2.10[ ), is the bulk free energy density. The 
latter exists only for t > 0. Thus, as a shortcoming of the 
Gaussian model, fex{t, L) can be defined only for t > Q 
although f{t,L), (2.9), exists for t < for the cases of 



antiperiodic, DD, and ND boundary conditions. 

The Casimir force -Fcas per component and per unit 
area divided by /cbT is related to /ox by 

Fcas(t,i)/(fcBr) = -^i^^^^^^. (5.2) 

For the subclass of isotropic systems, the asymptotic scal- 
ing form of its singular contribution is 



Fc^4t,L)/ikBT)^L-''X{i) 



(5.3) 



where, for i > 0, the universal scaling function X{x) 
is determined by the universal scaling function G{x) of 
the finite-size contribution to the free energy defined by 
(2.22) according to 



X{i) = {d- l)g{x) - i^-^i^^j^. (5.4) 

ax 

The surface contributions to the free energy density do 
not contribute to X{x). As an important consequence, 
finite-size scaling is found to be valid for the Casimir force 
for all b.c. in 1 < d < 4 dimensions, i.e., no scaling viola- 
tions exist for the Casimir force in the three-dimensional 
Gaussian model with NN, DD, and ND b.c, in contrast 
to the free energy density itself. 

As a consequence of ( 2.55 ) , the asymptotic scaling form 



of the singular part of the Casimir force becomes nonuni- 
versal in t he case of the anisotropic couplings (2.51). 
Then (5.3) is replaced by 



i^Cas,s(t,i)/(fcBr) 

= L-'^(Jx/J||)(^-i)/2x(t(L/eoa)'/'^), (5.5) 

where ^o,J_ is the amplitude of the correlation length 
(2.53b). Thus the Casimir force depends explicitly on 



the ratio of the microscopic couplings J± and Jy for all 
b.c, in agreement with earlier results for periodic [361138] 
and antiperiodic [41^ b.c. In the following we primarily 

discuss the isotropic case. 

For i > follow from with ( |420| and ( [XM] ) 



;f(NN)(~) ^ x(°D)(i) = 2-'^X(P)(4i), 

X(N°)(i) = 2-''x(^)(4i), 



(5.6a) 
(5.6b) 



and with (4.5b I 



X(")(i) = 2i-^X(P)(4i) - X^P\i). (5.7) 



dimensions by applying (5.4) to (4.5a). In three dimen- 
sions we utilize ( |4.9a ), its derivative 



Thus we only need X^'^\ which we obtain in 1 < d < 4 
by app] 
.hze (|4l 



dx 



47r 



In 1 



(5.8) 



and (5.7) to obtain the d = 3 scaling functions for peri- 



odic and antiperiodic b.c. for x > 0, 



X(P)(i) 



1 

TT 

X 

2^ 

_m 

TT 

1 



Li3(e-^) + ^Li2(e- 
In f 1 - e"^ 



-\/x\ 



X 

An 



Utt 
\n(l + e-^ 



(5.9a) 



TT 

X 

3C(3) i3/2 



Air 



12n 



+ 0(i'). 



(5.9b) 



The scaling functions for the other b.c. follow by employ- 
ing ([5^. 



At i = the critical Casimir amplitude [10] 

A = {d-1)-^X{0) (5.10) 
is obtained for 1 < d < 4 as 
A(P' - 2'^A(NN) - 2'^A(DD) ^ -^-''/2r(|)C(d), (5.11a) 



specializing for d = 3 to 



l)A(P), 



C(3) 
2tt 



^(a)^ 8^(ND) ^ 3C(3) ^ Q^^434g5^ 
87r 



(5.11b) 

-0.191313, (5.12a) 
(5.12b) 



The results (5.11) are identical to the results (5.6) and 
(5.7) of [Sj after setting = 1 (there is a misprint con- 
cerning the sign of Aper in (5.7) of [5]). The results for 
A*^P^ are also in agreement with Eq. (3.42) of [35] , 



As a consequence of (5.5), the Casimir amplitude 



Aaniso of thc anisotropic system (with Jj^ =^ J||) is 
nonuniversal and is related to A of the isotropic system 
(with J = J± = J||) for all b.c. by 

Aaniso = (Ji/J||)('^-i'/2A = (eo,±/eo,||)'-'A, (5.13) 

in agreement with earlier results for periodic [35] and 
antiperiodic [41] b.c. For d = 2, the right hand side 
of (5.13) is of the same form as found in [20] for the 
anisotropic Ising model on a two-dimensional strip with 
free b.c. 
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VI. CASIMIR FORCE IN d = 2 DIMENSIONS 

Our exact results for Xcauss in d = 2 dimensions are of 
particular interest in view of results for the exact Casimir 
force scaling functions Xising for an Ising model with 
isotropic couplings on a two-dimensional strip of infinite 
length and finite width L for free b.c. |23| and for periodic 
and antiperiodic b.c. in the recent work by Rudnick et al. 
[M] . Moreover, there exist earlier results for the Casimir 
amplitude at Tc of the two-dimensional Ising model with 
free b.c. and anisotropic couplings by Indekeu et al. j20j . 
This calls for a comparison with the corresponding Gaus- 
sian model results Xcauss- 



A. Isotropic case 

For the isotropic case, the two-dimensional Gaussian 
model results for periodic, antiperiodic, and DD b.c. are 



obtained from Eqs. g^, ( |43a| , Q, ( [5i3a| , and (jSJj) 

for d = 2. They read 



X. 



(p) 

Gauss 



1 

2^ 



dz 



3/2 



K{z)- 



X. 



(a) 

Gauss 
.(DD 



{x) 



54'aL(4i)' 

Gauss J 



X. 



(p) 

Gauss 



(6.1a) 

(6.1b) 
(6.1c) 



with the scaling variable x — (L/^q)^^"^^ — (L/S,)^, v — 
1/2, i>0. Above Tc, the corresponding Ising model re- 
sults for periodic, antiperiodic, and free b.c. read [231124] 



criticality) free b.c. in the Ising model correspond to DD 
b.c. in the Gaussian model. We see some similarity on a 
qualitative level: both the Gaussian and the Ising scaling 
functions are negative for periodic and for DD (or free) 
b.c, thus implying an attractive Casimir force whereas 
for antiperiodic b.c. they are positive implying a repul- 
sive Casimir force. On a quantitative level, however, the 
Casimir amplitudes at differ significantly, namely by a 



factor of two according to the exact results A'Q'g'uss(O) ~ 
2^1(0) = "-/6, 2A(fi_(G) = X(:4(G) = TT/e, and 
^Gausl(O) = 2X|1^^)(0) = -7r/24, obtained from (lIlOl, 



(5.11), and (6.2). Furthermore we note that, according 



to the dashed line in Fig. |4](a), Aq^^'^^^ has a weak max- 
imum above Tc, and correspondingly A^^^^g has a weak 
minimum above Tc, in agreement with Fig. 2 of |24j and 

Fig. 15 of tg. 

These results can be interpreted in terms of the two- 
dimensional Lp^ model which should be in the same uni- 
versality class as the two-dimensional Ising model. In all 
cases, the scaling functions at of the Gaussian model 
differ by a factor of two from the scaling functions at Tc 
of the two-dimensional (/S^ model. This indicates that a 
low-order perturbation approach in the fwo-dimensional 
(p^ model (in terms of a perturbation expansion with re- 
spect to the four-point coupling) is inappropriate, in con- 
trast to the situation in three dimensions to be discussed 
in Sec. |VII[ This is quite plausible since the fixed-point 
value of the renormalized four-point coupling in two di- 
mensions is quite large, i.e., far from the vanishing Gaus- 
sian fixed-point value. This is in line with the known 
fact that non-Gaussian fluctuations are generally larger 
in two than in three dimensions, as seen, e.g., from the 
bulk critical exponents. 



X- 



(a) 
Ising 



A, 



(ft-cc) 



(^l) 



(^l) 



1 

2^ 
1 

2^ 
1 

2^ 



duj q [tanh((7/2) - 1] , 
duj q [coth((7/2) — 1] , 



duj q 



[q+xi) e'J — {q~xi) e 
{q+xi)ei + {q-xi)e-i 



(6.2a) 

(6.2b) 

-1 
(6.2c) 



with q = \/x\ + uj"^. Here our Ising scaling variable xi = 

{L/S,o,if'''t = Lji with 1/ = 1 and ^o,/ = (8/3cJ)"^ [52] 
is related to the scaling variables A and x used in [53] 
and [23], respectively, by x\ = 2A = 2x (compare, e.g., 
with the isotropic limit of the correlation-length results 



in Appendix A 2 of [20]; see also Sec. VIB below). As an 
unexpected result, we find the surprising identities 



41ss((^/0') 



^Ishig(^/0' 



(6.3a) 
(6.3b) 



whose derivation will be presented elsewhere [44] . 

For a comparison of these results see Fig. [4] In 
Fig. |4j^b) , we have assumed that (asymptotically close to 



B. Anisotropic case 

We have extended the analysis of the isotropic Ising 
model by Rudnick et al. [21] for periodic and antiperiodic 
b.c. to the anisotropic Ising model on a square lattice 
with nearest-neighbor couplings Jy > and J_l > 0. 
This corresponds to the "rectangular lattice" of Indekeu 
et al. [2D] with the identifications of the couplings — 
2/3 J||, Ki = 2(5 J and ATs = 0. The corresponding 
bulk-correlation-length amplitudes above Tc follow from 
Appendix A 2 of [20] (for a = F) as [52] 



f -Ub T ^^^^ 

^''^ - 4 [^^-^^^ + smh(4/3cJ^) 

_ 1 / /3cJ|| 

^°-" = ir'^^ + sinh(4/3cJ„) 



(6.4a) 
(6.4b) 



with the ratio ( 2.60[ ), where we have used the condition 
sinh(4/3c J|| ) sinli(4/3c Jj_) = 1 for d = 2 bulk criticality. 
From Sec. Ill CI we obtain the nonuniversal Casimir force 
scaling function of the anisotropic Gaussian model for 
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FIG. 4: (Color online) Casimir force scaling functions (6.11 
and (6.2 1 of the Gaussian and Ising models in d = 2 dimen- 

to ' ' 



sions. (a) Scaling functions Xq 
(solid line) according 



■^cluaai^/O (dashed line) according to (6.3b I. (b) Scaling 



(6.3a I and Xj^{{L/0') 



Gauss 

functions X, 
line) . 



(DD) 



((L/C)^) (solid line) and xJ^(L/0 (dashed 



d> 2 above Tc 



^a„iso(i(i/eO,±)'^";^||,^±) 



(6.5) 



We have verified that this relation holds also for o? = 2 
dimensions (i) for the Gaussian model with periodic, an- 
tiperiodic, DD, NN, and ND b.c. and (ii) for the Ising- 
model for periodic and antiperiodic b.c. There is little 
doubt that it also holds for the two-dimensional Ising 
model with corresponding other b.c. Right at Tc this 
was established already in [501 for the case of free b.c, as 
noted in Sec. IlllC . 



VII. tp^ FIELD THEORY AT d = 3 

In this section we present the first (one-loop) step 
within the framework of the minimally renormalized cp'^ 
field theory at fixed dimensions 2 < c? < 4 [33l, ISGl for the 
calculation of the Casimir force scaling function in film 
geometry in the regime T > Tc for all five b.c. defined 

ITCl our 



m Sec. |IIj As noted in Sec. II C[ our Gaussian results 
for the various scaling functions G{x) and X(x) can be 
incorporated in such a theory based on the isotropic </3^ 
Hamiltonian (2.41 1. We emphasize, however, that we do 



not set uq equal to zero from the outset and that our one- 
loop treatment goes beyond the simple Gaussian model 
in that it includes the effect of the renormalized four- 
point coupling u via the exact exponent function Cr{u) 
(see Eq. (7.4) below), whose fixed-point value determines 
the exact (non-Gaussian) critical exponent i'. The one- 
loop approximation manifests itself only in neglecting the 
(two- loop) 0{u*) contribution to the amplitude function 
of the free energy density. Such a treatment has recently 
been presented in Sec. XA of [36| for the case of cubic 
geometry with periodic b.c. For the specific heat in film 
geometry with Dirichlet b.c. a corresponding treatment 
was given in [T2]. As suggested by the earlier successes 
[T2l [TTI [36] , the minimally renormalized (p"^ theory at 
fixed d is expected to constitute an important alternative 
in the determination of the Casimir force scaling func- 
tion in comparison to the earlier e expansion approach 
[SI m [15]. Our quantitative results to be presented in 
Fig. [5]below will support this expectation. Other fixed-d 
renormalization schemes are, of course, conceivable which 
would lead to the same one-loop results at d = 3 as ob- 
tained in our approach. We believe, however, that the 
fixed-d minimal subtraction scheme has considerable ad- 
vantages in extending the finite-size theory to two-loop 
order and to the temperature regime below Tc [HI [36] . 

As a temperature variable we use the shifted parame- 
ter ro — tqc = aot, where tqc = — 4(n -I- 2)uo J^'*'' fe~^ is 
the critical value of tq up to O{uo). In the following we 
sketch the relevant steps of calculating the singular part 
of the minimally renormalized free energy density in one- 
loop order for film geometry with periodic or antiperiodic 
b.c. After subtracting the regular bulk part up to lin- 
ear order in tq — tqc and performing the limit A — > oo 
at fixed To — tqc we obtain the bare one-loop expression 
of the remaining part 6f of the free energy density per 
component divided by ksT in 2 < d < 4 dimensions as 

<5/(fo - rocUo, L) = (ro - tqc)''^^ 

+ rt((ro-roc)i')+OK) (7.1) 
for periodic and anti perio dic b.c^ where G^^'^iy) and 



Q'^'^\y) are given by (4.5a I and (4.5b). The renormal 
ized parameters r and u are defined in the standard way 
[36] as r = Z^^{rQ — roc) and u = fi^^ A^Z:^^ Z^uq with 
an inverse reference length /i = CcT^j where is the 
correlation- length amplitude above Tc- The additively 
renormalized counterpart /r of Sf is defined as [36] 



fR{r,u,L,fx)^Sf{Zrr,^fZuZ^ u, L) 



1 



fi 'r^AdA{u,e), 



(7.2) 



where A{u,e) = —2n/e + 0{u) is the additive renormal- 
ization constant of the minimal renormalization scheme. 
After integration of the renomalization-group equation 
(see Eqs. (5.6) and (5.7) of [35]) and with the choice 
^2/2 _ ^^j^ q£ ^j^g gQ.^ parameter Z, the finite-size part of 
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fn is then given by 

/fl(r, u, L, ii) - Mr, u, oo, ^^) = L-''g{r{l)L') + 0(u(l)), 

(7.3) 

with the effective temperature variable |33) 



r(0=rexp( f CrW))'^ 



(7.4) 



This variable contains the field-theoretic function Cr(w) — 
IJ,{dfi,\n Z~^)o 03] with the fixed point value Cr(w*) = 
2 — i'~^. Asymptotically (^ — ?> 0), this leads to the scaling 
form of the singular part of the excess free energy per 
component in one-loop order 

/e.,s(i, L) = L-^GiL^/e) + 0{u*) (7.5) 

for periodic and anti perio dic b .c, with t^(p)(y) and 
G^^Hy) given by (4.5a) and (4.5b) where now = ^o^""^ 
is the correlation length above with the exact critical 
exponent v of the (d, n) universality class. 

A complication appears to arise at d = 3 for NN and 
DD b.c. because in these cases the limit A — > oo at 
fixed rg — tqc does not exist which in the dimension- 
ally regularized form of the free energy density shows up 
as pole terms at d = 3. These divergent contributions, 
however, are restricted only to the surface part L~^fsf 
which does not contribute to the Casimir force. Thus 
fcx{t, L) — L^^fsi{t) is well behaved at d = 3 in the limit 
A — >■ oo at fixed rg — tqc- Consequently there exists no 
problem in the calculation of the finite-size part scaling 
function Q and of the Casimir force scaling function X in 
the framework of the minimally renormalized theory at 
fixed d = 3 for NN and DD b.c. This holds also for ND 
b.c. The resulting asymptotic scaling form in one-loop 
order is 

/e,,,(i, L) - L~^fsiAt) = L-^GiLye) + 0{u*), (7.6) 
for NN, DD, and ND b.c, where G^^^Uv). and 



G^^^Hy) are given by (4.20) with (4.5a), and by (4.36) 



with (|4.5b[ ), respectively. Again, the correlation length 
^ in ( |7.6[ ) contains the exact critical exponent i> of the 
(d, n) universality class. The results (7.5) and (7.6) can 
be applied directly to d = 3 dimensions. T he C asimir 
force scaling functions A(L^/(^^) follow from (5.4). Thus 



we are in the position to perform a reasonable comparison 
both with MC data for the three-dimensional Ising model 
[B] and with earlier e = 4 — d expansion results [S', TF of 
the (p'^ theory evaluated at £ = 1 as well as with the very 
recent result of an improved d = 3 perturbation theory 
[T7] (in a L| X L slab geometry with a finite aspect ratio 
p = L/L\\ — 1/4) for periodic b.c. This comparison is 
one of the central results of this paper. 

The comparison is shown in Figs.[5]ja)-(j) as a function 
of the variable L/<^. The solid lines represent our one-loop 
results. The e expansion results are represented by the 
dashed lines (two- loop e expansion [Q) and by the dotted 



lines (improved e expansion [H]), respectively, and the 
result of the improved d = 3 perturbation theory [17] 
is represented by dot-dashed fines in Figs.j5][a),(b). In 
the large-L/^ regime (not shown in Fig. [sF, the solid 
and dashed lines have an exponential approach to zero 
and differ very little from each other in all cases. This 
statement holds also for the dotted and dot-dashed lines 
for the case of periodic b.c. (Figs. [5] (a), (b)) but not for 
the dotted lines for the case of NN b.c. (Figs.[5](e),(f)), 
where the e expansion result of [15j breaks down in the 
large-L/^ regime. Also shown are recent MC data [B] 
for the three-dimensional Ising model (in a L| x L slab 
geometry with the aspect ratio p = l^jl^w =1/6 and with 
L = 20) for the cases of periodic and DD b.c. in Figs. [5] 
(a),(b),(g),(h), respectively. 

For periodic b.c, our one-loop result for G^'^^L'^ / S,'^) 
and A(p)(LV^2) (goij^ y^^^^ Figs. [H]; a) and (b)) is in 
remarkable agreement with the improved e expansion re- 
sult of [15j (dotted lines in Figs.[5]ja) and (b)) at (i.e., 
i/^ = 0) and in the large-L/^ regime. The slope of our 
one-loop result at is in better agreement with the slope 
of the MC data than the slopes of the e expansion results 
at Tc. In particular, there is no artifact of the one- loop 
result of the fixed d = 3 theory such as the minimum of 
the two-loop £ expansion result above Tc shown by the 
dashed line in Figs. [5ja) and (b). This suggests that, 
for periodic b.c, the d = 3 approach is a better start- 
ing point of perturbation theory than the e expansion 
around d — A dimensions. This is consistent with recent 
findings of finite-size effects in cubic geometry (Fig. 5 
in 36J) and in finite-slab geometry (Fig. 4 in [IT]). For 
antiperiodic b.c. (see Figs.[5jc) and (d)), there is surpris- 
ingly good agreement between our one-loop d = 3 result 
for G'^'^KL^/i'^) and X('^)(LVC^) and the two-loop e ex- 
pansion result [Hj. This is quite remarkable in view of 
the fact that the computational effort in obtaining d = 3 
one-loop RG results is considerably smaller than that for 
deriving two-loop e expansion results. 

For NN, DD, and ND b.c, there are considerable dif- 
ferences between our one-loop results of the fixed d = 3 
ip'^ theory and the two-loop e expansion results as shown 
in Figs. [5]je)-(j). For NN b.c, however, the improved e 
expansion result of |15| at Tc is not far from our one-loop 
result, but for L/^ > our d = 3 result does not agree 
with the strong increase of the e expansion result of [15j . 

In summary, the fixed d = 3 theory yields reasonable 
results already in one-loop order. It would be a reward- 
ing task to perform a two-loop calculation of the fixed d 
theory for all b.c. as well as to proceed to higher-orders 
of the £ expansion for DD and ND b.c. Also MC data 
for the cases of antiperiodic, NN, and ND b.c. are highly 
desirable for a comparison with the predictions shown in 
Fig. [5] in order to clarify the reliability of the different 
perturbative approaches. 



23 




x(p) 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 



0.15 



0.10 



0.05 



0.00 L 



(c) 



1 



g(NN) 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 




0.00 



-0.01 



DD) 



-0.02 



0.020 r 




0.00 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 

m 



0.015 
^(ND) -.^^ 
0.010 



0.005 



(i) 



0.000 L 



0.0 



0.5 



1.0 



1.5 



2.0 



2.5 



3.0 




FIG. 5: (Color online) Scaling functions Q and X of the finite-size part of the free energy density and of the Casimir force, respectively, 
for T > Tc in a three-dimensional film of thickness L with isotropic interactions fo r the variou s b.c. a s a fun ction of th e scaling argument 
L/^. Solid lines: one- loop results of the fixed d = 3 ip'* theory according to Eqs. | |4.26| , | |4.46| , ^T^ , and for and 

according to Eqs. i 5.6 I, ( 5.9 i, and ( 5.4 i for X((L/^)^). Dashed lines in left panels: two-loop e expansion results S+pcrli/^); 0+apGr(i/5), 
0-I-SB,Sb(^/5)i 0+0,0(^7?)) and B_|_o SB(i/5) at e = 1 for n = 1 for periodic, antiperiodic, NN, DD, and ND b.c, respectively, according 
to Eqs. (6.12) and (6.13) in [5]. Dashed lines in right panels: two-loop e expansion results i9+per(i/^), i?+aper(i'/5), i?+SB,SB (i'/?)i 
i?+0 o(^/?)i and i?+o Sb(^/?) at e = 1 for n = 1, obtained through Eq. (3.9) in 9^. Dotted line in (a): improved e expansion result 
e(P<='')(L/5) at £ = 1 for n = 1 according to Eq. (4.57) in [l5l, compare Fig. 6 of [15]. Dot-dashed line in (a): F'=''((L/5)l/'', p = 1/4) in 
fixed d = 3 according to Eq. (17) in [ij. Data point in (a): From MC result i?p(0) = -0.3040(4) in [6]. Dotted hne in (b): improved e 
expansion results (L/^) for n = 1, obtained through Eq. (1.7) in [15]. Dot-dashed line in (b): X((L/g)l/"', p = 1/4) in fixed d = 3 

according to Eq. (19) in 1171 . The L/^ = data point in (b) is twice the L/g = data point displayed in (a). Other data points in 
(b): MC results from Fig. 15 in [6] for L = 20 and p = 1/6. Dotted line in (e): improved e expansion result e(=P'=P'(L/5) at e = 1 for 
n = 1 according to Eqs. (4.47), (4.50), (4.53), and (4.56) in I15| . compare Fig. 8 of |15 | (where the variable on the abscissa should read 
\/L/^oo = x/L)- Dotted line in (f): improved e expansion results E:(=P''=p)(L/0 at e = 1 for n = 1, obtained through Eq. (1.7) in [TSl . 
Data point in (g): MC result for L = 20 and p = 1/6 by dividing the L/^ = 0-result displayed in (h) by 2. Data points in (h): MC results 
for L = 20 and p = 1/6 from the inset of Fig. 13 in [6] (where the e expansion line is misrepresented). 
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VIII. DIMENSIONAL CROSSOVER: 
HEAT 



SPECIFIC 



In the following, we present an explicit study of the 
dimensional crossover in the Gaussian model from the 
finite-size critical behavior near the d-dimensional bulk 
transition at to the (d— l)-dimensional critical behav- 
ior near the film transition at the critical temperature 
7c, film (-^) < Tc of the film of finite thickness L. (The 
equality sign holds only for periodic and NN b.c, whereas 
Tc,fiim(^) < Tc for antiperiodic, DD, and ND b.c.) The 
most interesting candidate for this study is the divergent 
specific heat C(t,L), whose d-dependent critical expo- 
nent a ch anges from abuik = (4 — (i)/2 near bulk Tc (see 
Eq. ( pSp to amm = [4 - (d - l)]/2 near Tc3im(i) m 
d < 4 dimensions. We shall verify that the film critical 
behavior of a d-dimensional system corresponds to that 
of a bulk system in d — 1 dimensions for all b.c. except 
for antiperiodic b.c, where an unexpected factor of two 
appears due to a two-fold degeneracy of the lowest mode 



as noted already in Sec. 



above. It remains to be 



seen which effect this feature may have in non-Gaussian 
models and in the (p'^ theory. 

The dimensional crossover behavior is particularly sim- 
ple in the Gaussian model because the correlation-length 
exponent = 1/2 is independent of the dimension d and 
the correlation-length amplitude is the same both for the 
bulk and the film critical point. To the best of our knowl- 
edge, this crossover behavior has not been presented in 
the literature so far. Nevertheless it is worthwhile to 
study the exact Gaussian crossover behavior for various 
b.c. as it provides part of the mathematical basis also 
for the crossover behavior in the more complicated mean 
spherical model in film geometry in 3 < d < 4 dimensions 
that we shall study in a separate paper [14] . 

According to 



specific heat reads 



2.9 1 and (2.251, the expression for the 



CXt,L) 



TrL 



T. 2-^ 

1 

E 



p 



+ Jp4-1 + Jq 



(8.1) 



For small t and large i/a, the specific heat can be de- 
composed into singular and nonsingular parts as 



C{t,L) = C,{t,L) + CUt.L). 



(8.2) 



The first term on the right hand side of (8.1 1 provides the 
leading singular contribution to Cg, whereas the second 
term yields only subleading corrections. The non-scaling 
structures of the type discussed in the preceding sections 
(for NN, DD, and ND b.c. in d = 3 dimensions) appear 
only in the subleading corrections, whereas the leading 
part of the first term of (8.1) is in full agreement with 
the finite-size scaling form (for the subclass of isotropic 
systems) 



(8.3) 



as noted already in [3T] for the case of free (DD) b.c. For 
the Gaussian model, the scaling structur e (|8.3[ ), together 
with the critical exponents ( |2.12 ) and ( |2.28[ ), holds in 
1 < d < 4 dimensions for all boundary conditions. If 
the finite-size scaling function T{x) of the free energy 
density, ( 2.18 1, exists, the universal scaling function C{x) 



is related to it by 



C{x) 



d^Tji) 
dx^ 



(8.4) 



For the simplest anisotropic case, discussed at the end of 
Sec. |II C[ the corresponding nonuniversal result can be 
derived from (2.55). C{x) exists also in d = 3 dimensions 
for the cases of NN, DD, and ND b.c, where F{x) does 
not exist for all x. 



From (2.201 follows the bulk singular part 



with the universal bulk amplitude 



YcAo^^^^/-, (8-5) 



Yc.d = - 



d{d-2)^__ r(4^) 



2(47r)'i/2 



(8.6) 



which is valid for d>0, d7^4;6,8,... dimensions, even 
though (2.20) does not hold for d = 2. This implies in 



two and three dimensions 



1 



8^^o ^ ' 



SO 



1 



a,s(t) = 

167r''" " " 167r 
The finite-size scaling functions read 



d: 
d: 



?.7a) 



c(p)(i) = -/cf (^) + yad^(''-4'/^ 

C(^)(i) = -K.f{x + TT^) + YcAi + 7^^)(''-'^/^ 



i > 0, 

X > — TT^, 

i > 0, 

X > — TT^, 

X > -(7r/2)2, 



(8.8a) 
(8.8b) 
(8.8c) 
(8.8d) 
(8.8e) 
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with 



A 



(N) _ 
C,sf ~ 



-A. 



(D) 
C,sf ■ 



_ (d-i)(d-3) ^^(N) _ r(^) 



8(47r)('*-i)/2' 



(8.9) 



as follows from (18^,(|41|, ( |4.16| ), and ( |4.40| ). The func- 
tions lCd{y) are listed in Appendix |Dj They are given 
by the second derivatives of the functions Td defined in 
Sec. 



IV 



I.e. 



Eqs. (8.8) are valid in 



1 < d < 4 dimensions. Eq. (8.9) agrees with Eqs. (125) 
and (126) of [3T] for the case of DD b.c. (apart from a 
factor two due to a different definition of ^c^sf)- 

The scaling functions (8.8 1 are valid not only near 



X = 0, but also near the film transition at Tc^fiini(i) < Tc 
for antiperiodic, DD, and ND b.c, i.e., near x — — tt^ 
or i = — 7r^/4, respectively. This means that they 
provide an exact description of the crossover from the 
d-dimensional to the (d— l)-di mensiq nal critical behavior 
of the specific heat. The result (8.8d) for C'^^'(x) agrees 
with the result for C(y, 0) of Eqs^T24)-(126) in with 
the identification x — . 



Similar to (3.6), at finite L, we define the film specific 
heat Cfiim (heat capacity per unit area divided by fee) as 



Cfilm(ro,L) =iC(t,L). 



(8.10) 



One expects that, asymptotically (^fiim ^ L), the 
film critical behavior of a d-dimensional system cor- 
responds to that of a bulk system in d — 1 dimen- 
sions. We indeed obtain from (8.8) and (8.31 for small 
['"0 — rQcfLimiL)]/ J <^ L^^ the singular part of the film 
specific heat for finite L for all b.c. 

Cfilm,s(^film, i) = 

ic\d-i ,(rf-5)/2 _ Ycd-i periodic, NN 



Co 



d-1 ''film 



?film ' 



^4 --mm ' DD, ND b.c, 

(8.11a) 

21c.d-i (d-5)/2 _ 2Yc,d-i r,_d antineriodic b c 

^d~l film ^ c4 ^film ' dnupcllOUlL U.C, 

(8.11b) 



with tfiim = t — tc.fiim{L), in agreement with the bulk 
critical behavior in d—1 dimensions (compare (8.5) and 
(8.7), as expected on the basis of universality. An ex- 



ception is the additional factor of 2 for antiperiodic b.c. 
which is a consequence of the two-fold degeneracy of the 
lowest mode, as already noted in Sec. Ill (see Eq. (3.8 1). 

Since the crossover behavior is qualitatively similar in 
all dimensions 1 < d < 4 we confine ourselves to illus- 
trating only the example of DD b.c. in d = 3 dimensions 
which is obtained from (8.8d) as 



1 / coth-\/i 1 
167r I ^/^ X 



X > —IT 



i.l2) 



(Eq. (|8.12 ) fo llow s also from 
and (4.261 or (E2|.) Its asymptotic behavior is 



.1) together with (4.29) 
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FIG. 6: (Color online) Solid line: Double- log arithm ic plot 
of the specific-heat scaling function C''^'^'(x), (8.121, of the 
Gaussian model with DD b.c. in three dimensions as a func- 
tion of i -I- TT^, with X = t{L/^of, t={T~ Tc)/Tc. Dotted: 



asymptotic behavior at small x + tt^ > according to 1 8. 13a I 
and (8. 14a I, displaying the divergent critical behavior near 



Tc,fiim(i) < Tc with an exponent a — 1. Dashed: asymptotic 
behavior at large x according to (8.13c I and (8.14c I, display- 



ing the three-dimensional bulk critical behavior above bulk Tc 
with an exponent a — 1/2. Near a: = corresponding to bulk 
is an analytic function with the finite amplitude 
C{DD)(o) = l/(487r). 



C(DD)(5 



1 



+ 0{{x + nT), 



8t:{x + 7r2) 



1 

48^ 
1 



7207r 



1 



0{i 

^ +0(e-^/^). 



< i 4- TT^ < 1, 
(8.13a) 
|i| < 1, (8.13b) 

5>1, (8.13c) 



167rV5 167rx 
which implies the corresponding asymptotic behavior 



1 

1 



f-2r-l,-l 
?0 ^ %lm' 



487r- 

J_^-3,-l/2 



< i/Cmm « 1, (8.14a) 
T^Tc, L/d>l, (8.14b) 
L/^ > 1. (8.14c) 



This indeed represents a two-dimensional critical behav- 
ior near Tc^fiim with the exponent a = 1 (compare (8.7a)), 
a three-dimensional finite-size critical behavior at Tc with 
I' = 1/2, a = 1/2 (compare Eq. (8.3)), and a three- 



dimensional bulk critical behavior above Tc with the ex- 
ponent a = 1/2 (compare Eq . (|8.7b )), respectively. The 
crossover is illustrated in Fig. |6J Similar illustrations can 
be given for the other b.c. 
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Appendix A: Film Thiclcness 



Here we derive the expression ( 2.48 1 for the thickness L 



of the transformed isotropic film. This thickness is given 
by L times the height of the cZ-dimensional parallelepiped 
spanned by the transforms x'^ = X^^^^Uxi of the orthog- 
onal unit vectors Xi^ i = 1, . . . ,d over the surface given 
by the (d— l)-dimensional parallelepiped spanned by x'^, 
i = 1, . . . , d—1. Since the volume of a parallelepiped is 
given by the product of one of its surface areas and the 
corresponding height, we may write 



L = {Vd/Vd-i)L, 



where the volume of the d-dimensional parallelepiped 
has been denoted by Vd and the surface area is given 
by the volume Vd-i of the (rf— l)-dimensional paral- 
lelepiped. The TO-dimensional volume Vm of a par- 
allelepiped spanned by the m vectors Vi,. . . ,v„i in an 
n-dimensional space with n > m is given by Vm = 

1/2 

(det V'^V) , where V is the n x m matrix whose 
columns are the vectors vi,. . . ,Vm- Thus we obtain 

Vd = [det(A-i/2(7)^(A-i/2(7) 

d 

'1/2 4-ni/2 



1/2 



= (dot 



-ni/2 



i=l 

where A^^ may also be written as 



(A2) 



(A3) 



V 



and 



Vd-i = (det[A-i/2(7]^[A-i/2j7] 
= (det[A-i/2]T[^-i/ 



1/2 



{det[[A-']]f\ (A4) 



where [A~^/^C/] and are the d x (d— 1) matrices 

that result from removing the last column from the ma- 
trices X~^/^U and A^^^'^, respectively, and where [[A^-'^]] 
is the (rf— 1) X (d— 1) left upper part of A^^. Combining 
(|Al|), (|A2|), and (tAil) gives ^A8\. 



For the specia 



geometry was investigated, of which the film geometry is 
only a limiting case; (ii) we perform an exact separation 
of the surface contributions at an early stage of the cal- 
culation and reduce the remaining computations for all 
b.c. to the case of periodic b.c. 

First we cons ider the isotropic case. We start from 

J,°°dyy~He'y 



(2.9) and (2.10) and use \nz 
to write the excess free energy as 



-yz 



/ex(i,i) 



1 



2a' 



dy 



e-y^o/^Biyf-^ABr^iy), (Bl) 



(Al) with fo defined after Eq. (2.31) 



ABN{y)^Biy)~BNiy), 



B{y) 



1 

2ti 



dip exp[— y(l — costp)], 



(B2) 



(B3) 



SAr(y) = -^^exp[-y(l-cos9™a)], (B4) 
where the sum ^ runs over the wave numbers given 



in (2.4). The quantity B{y) = limAr_>.oo BMiy) in (B3) is 
identical to B{y) from (2.33). By rearranging the sums 



it is possible to express B^^y), B^^\y), B^^\y), and 



b'n^^v) terms of B''^'(y). For example, for DD b.c. 



?(p) 



B^Y'-'iy) 




exp 



—y 1— cos 



7r(m+l) 
iV+1 



2(JV+1) 



i)(y)' 



(B5) 



where in the first step we have exploited the symmetry of 
the cosine about tt, while in the second step m = — 1 and 
m — N terms have been added to and subtracted from 
the sum and subsequently m has been renamed m — \. 
Similar rearrangements can be performed for the other 
nonperiodic b.c. and we obtain the exact relations 



case where A is diagonal, we obtain AB^iy) = 2AS^^(i/) - AB]^\y) 



Vd = ntiAJ^'^' and Vd-i - nti A^'^' and thus 



Vd/Vd-i - A 



-1/2 
dd 



Appendix B: Free energy 



AB 



AB 



(NN) 
N 

(DD) 
N 



Here we derive fs{t, L) of the Gaussian lattice model in 
1 < d < 4 dimensions for film geometry with the various 
b.c. for both the isotropic case and the anisotropic case. 

While we follow in spirit the derivation given for DD 
b.c. in [3T], two simplifications arise: (i) In pi] a slab 



AB 



(ND) 
N 



(y) 
(y) 



(y) 



l+e-^y-2B{y) 



2N 



1+- I AB 



(B6a) 
(B6b) 

2('Ar+l)(2^)' 

(B6c) 



e-^y - B{y) 



2N 



2N 



AB. 



(B6d) 



This leads to the exact representation 



27 



/cx(i,L) 



2/sf(0 



1 

2^ 



'N 

4w( 



(1 
(1 



1 ^ AB, 



2(Af+i)(y)' 
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periodic b.c, (B7a) 

antiperiodic b.c, (B7b) 

NN b.c, (B7c) 

DD b.c, (B7d) 

ND b.c, (B7e) 



for arbitrary L = Na with AB^^ from (|B6a|). For NN, 



DD, and ND b.c, the surface contribution 2fsi{t)/ L origi - 
nates from t he firs t term on the right hand sides of ( B6b I , 

and (B6d), respectively. 



(|B6c 



It is g iven by / f (Q, 

f^^'{t), and f^"'^\t) provided in ( |4.1l[ ) and (2.14|. 
The remaining tasks are (i) to determine the large- A'' 



Extending the lower integration limit in /2 to leads only 
to exponentially small corrections. Changing the integra- 
tion variable accord ing t o z — 2y(7r/iV)^ gives for peri- 



odic b.c^the resuh (2.221 with g^P\x) in (4.5a 



Due to 



(|B7bl)-([B7e| , we confirm (|2.22|) also for the other b.c. un 



der consideration here (the surface terms are absent also 
for antiperiodic b.c), with the g{S:) provided in (4.5b 
(4.20), and (4.361, with specializations to d = 3 in (4.9 



(|4^, and (|446|). 

Now add to L~'^Q{x) the bulk singular part of the free 
energy /b,s from ( 2.20 ) and, for NN or DD b.c. and d 7^ 3, 



the surface singular part from (4.12 ) (the correspond- 



ing part for ND b.c. vanishes, see Sec. IV D). Observing 



( 2.18 ) with Ci given after ( 2.19 ) leads to the scaling func- 
tions (|4T]), ( |4.16[ ), and ( |4.4G| ). 

For d = 3 and NN or DD b.c, we add to L-'^Q{x) with 



behavior of B^^\y) for periodic b.c. and (ii) to translate singular parF^rom ( |4.23[ ) to obtain tliel-esuhs ( |4.25 l 
the result to the other b.c. For the first task it is useful 
to distin guish the regimes < y yo and y ^ j/o in the 
integral (Bl|, with j/o chosen such that 1 <C yo ^ 
Thus, for periodic b.c, we separate 



Gix) from (4.261 the bulk singular part of the free energy 
/b s from (2.20) with from after (2.34) and the surface 



Inserting the small- 2: expansion of K(z) into (4.5a 
with ( |4.2a ) gives 



/ex(t,i) = ^(/l+/2), 



(B8) 



dz 



2Z /7r\ 

TT \ Z J 



(d+2)/2 



n=l 



(B13) 



r I dy _ 

/i,2 = / — e 

'1,2 y 



■yo/2 1 



Biyr-'AB)!^'{y) 



(p), 



(B9) 



with J-^ = and = , corresponding to (AlO) 
and (All) of 31,. As shown in (A12)-(A17) of [S], the 
large- iV dependence of AB'j^\y) in the regime < y 
yo is of 0(e~^), thus /i yields only exponentially small 
contributions. 

Now consider y ^ yo with yo ^ 1. Rewrite the 
sum over qm in (B4) for periodic b.c. by letting m 
run over m = —N/2, . . . , A^/2 — 1 for even N and 
m = -(7V-l)/2,...,(iV-l)/2 for odd N. Then only 



l^mol ^ 1 can lead to contributions to B^^\y) in (|B4l) 



that are not exponentially small and we may expanc 

exp[-y(l - cosq^a)] 

= exp {-iy(g„a)2 [1 + 0(((7„d)2)] } . (BIO) 

Correspondingly, we obtain in the regime y ^ yo , apart 
from exponentially small corrections, 



1 

N 



+ CX3 

E 



exp 



(Bll) 



For the evaluation of /2, we use (Bll) and keep only 
the first term of 

B{y) = -^[l + 0{y-')]. (B12) 



For large i, the right hand side is dominated by the first 
term of the sum. The remaining integral may be eval- 
uated in a saddle point approximation. This yields, to- 
gether with (4.5b|, (4.20), and ( 4.36 1, t he ex ponent ial de - 
cay of the non-surface terms in (14. 71), (4.22), and (4.38). 



For the anisotropic case, we consider two different cou- 



plings J± and Jii, see (2.51). Then (Bl) is replaced by 



Ta'^ Jo y 



~, , e-y'"'-/'B{i.mj^)yf-'ABN{y), (B14) 



with fo,_L = f'QcP' / {2J±). For the leading singular finite- 
size terms only the leading large-y behavior of B{y) given 
in (B12l matters, see the derivation of (4.5a) for periodic 



b.c. above and its translation through (B7b)-( |B7e ) to the 
other b.c, manifested in (|4.5b|, (|4.20i, and (|4.36|). The 



same is true for the leading singular surface terms as may 
be inferred from the derivation of (4.12) for d 7^ 3. Thus 



the factor J\\/ J± in the argument of B in (B14) leads to 



an additional factor ( Jj^/ Jy )('* ^^/^ in front of the leading 
singular contributions to the excess free energy /ox(^7 L). 
Eqs. ( |2.31| -( p32| ) are replaced by 



hit) = 



Wa{z,w) 



2d.'^ 



In 
dy 

y 



+ T^d(roa, J||/J±) 



(B15) 



-v/2 _ e-'^y/^B{y)B{wy) 



(B16) 
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Thus, because of the argument {J\\/J±)y of B, the same 
factor {J±/ J\\Y'^~^^^'^ appears in front of the leading sin- 
gular bulk part. Since the temperature dependence en- 
ters only through the parameter fo,_L — ^^/di where ^± 
is the correlation length (2.53b), it is straightforward to 



confirm ( 2.55 1 for all b.c. on the basis of these properties. 



Appendix C: Comparison with Ref. [S] 

In Sec.^ we stated the identity oiQ{x) with the func- 
tions 6^'(j/+) used in [S] (as noted for periodic and an- 

for NN and DD b.c. af- 



tiperiodi c b.c . after Eq. (|4.5 
ter Eq. (|420]), and for ND 



b.c. after Eq. (4.36 1). For 



periodic b.c, this equivalence follows from using the ex- 
pansion (3.38) in [3S] in terms of Bessel functions, which 
provides another representation of With i = y^, 

we obtain 



1 

'271 



dz (7r/2)(''+i)/2e— 4/(2-)^ \k{z) - 



= --E / dz{'Klzfl^+^e-^-yM'^^^^\-^"^"l^ 

^ ^ 1 



rf/2 Kd/2{ny+) 



e 



(47r)('i-i)/2r( 

■(4^)(<i-i)/2r(rf±i) 7^ 
(1) 



dziyZ^ 



X')('i-i)/2g-«2y+ 



dz 



(z2-l)('^-l)/2 

e2a+ - 1 



Hper(Z/+), (CI) 

valid for 2/4. > 0. The other identities between Q{x) and 

0^1^ ''(2/+) follow similarly. They may also be derived by 
showing the identities (4.5b), (4.20), and (4.36) for the 

fimctions O^-* as represented in [5]. 

Appendix D: Functions 

For \z\ < 1, the poly logarithms are defined by Lii,{z) ~ 
z'' /k'^ and for \z\ > 1 by their analytic continua- 
tion. They are analytic in the complex plane except at 
z = I and except for a branch cut that we take along 
z oo[. We need Lii,{z) for = 1, 2, 3. Well known re- 
lations are Lii,(l) = C{iy) for Rei/ > 1, Li^(-l) = {2^-"- 
l)C,{u), Lii(z) = — ln(l — z), and Li[,(z) = Li^_i(z)/z for 
z ^ [l,oo[, where ({v) = YlV=i^~'' Riemann's zeta 
function. Combining them we may write for z > 

Li3(±e-") + zLi2(±e-") 

= i(l±7)C(3)+^'da;xln(lTe-^), (Dl) 
which is needed for App. [E| 



The various functions /C^ appearing in (8.8) read 

POO 



ir(DD)/ \ 



1 



327r3 



dz{7r/zf~'^^\-^y/<^'^^\ 







(D2a) 



K{z) — \J tt/z 

— r 

32^3 7o 

|e^/-* [A'(z/4) - K{z)\ - ^/^z'^ , (D2b) 



dz{^/z)'^^-^^'^e-^-y'^^^^\ 



1 



2^+17^3 



dzi^/zi'-^^'^e-^y/^^ 



\e''[K{z)-l]- ^/V/^z + lY (D2c) 
and/Cr^(y) ^ 2^-'^Kf {Ayl JCf'Xy) = 2^'^Kf{Ay). 

Appendix E: Analyticity properties 



Here we s how t hat the d = 3 expressions ( 4.8b[ ) for 
J-(=^)(i) and jliol for L^f^^^\t,L) + (87r)-ii ln(L/a) 
are analytic for x > —t t'^ an d are finite and real at a; = 
— 7r2. Combining (4.9b) and (Dl|, we obtain for x > 



- - 

47r 1 2 



1 

2^ 



rC(3) 



/ dzz\n{l 
Jo 



C(3) 



^3/2 



dx' In 



2cosh(Vx'/2) 



(El) 



Since cosh(vi') is analytic and pos itive at i' = 0, the 
analytic continuation of J^*^^^(x) in (4.8b) from positive 



X to other x is analytic at x = 0. Inspection of (El) 



shows that the in tegra l there is also analytic for all other 

^2 



> — 7r2. Thus (4.8b) is analytic in x for all x > — tt' 



Computing the integral in (El I for x = —tt'^ and combin 



ing the result with (4.8b) gives the result (4.10). Com- 
bining (4.26) and (Dl), we obtain for i > 



1 

16^ 
1 

16^ 



C(3) 



/ dzz\u{l-e-') 
Jo 



C(3)- TT5^/'+i(lni-l) 



2/%i'ln^^-h(^) 



(E2) 



Since sinh(vi')/vi' is analytic and positive at x' — 0, 
the analytic co ntin uation of (4.30) is analytic at x = 0. 
The in tegra l in ( E2 ) is also analytic for all other x > —tt"^ . 
Thus (|4.30t is analytic in x for all x > -n^. Using (|E2|) 



to ex pand (4.30) around x — — t t for x > — tt gives 
(4.31), with a finite value of (4.30) at a; = —tt'^. 
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